re ty 
hi ) i A a] 


The, . 
ilosophical Magazine 


FIRST PUBLISHED IN 1798 


A Journal of Theoretical 


Experimental and Applied Physics 


Vol. 3 April 1958 No. 28 
Eighth Series 


£1 5s. Od., plus postage 
Annual Subscription £13 10s. 0d., payable in advance 


Printed and Published by 


TAYLOR & FRANCIS LTD 
RED LION COURT, FLEET STREET, LONDON, E.C.4 


THE PHILOSOPHICAL MAGAZINE 


Editor 
Professor N. F. Mort, M.A., D.Sc., F.R.S. 


Editorial Board 


Sir LAWRENCE Braaa, O.B.E., M.C., M.A., D.Sc., F.R.S. 
Sir Grorar Tuomson, M.A., D.Sc., F.RB.S. 
Professor A. M. TynpDALL, C.B.E., D.Sc., F.R.S. 


AvutHors wishing to submit papers for publication in the Journal should 
send manuscripts directly to the Publishers. 


Manuscripts should be typed in double spacing on one side of quarto 
(8x10in.) paper, and authors are urged to aim at absolute clarity of 
meaning and an attractive presentation of their texts. 


References should be listed at the end in alphabetical order of authors 
and should be cited in the text in terms of author’s name and date. Dia- 
grams should normally be in Indian ink on white card, with lettering in 
soft pencil, the captions being typed on a separate sheet. 


A leaflet giving detailed instructions to authors on the preparation of papers 
is available on request from the Publishers. 


Authors are entitled to receive 25 offprints of a paper in the Journal free 
of charge, and additional offprints can be obtained from the Publishers. 


The Philosophical Magazine and its companion journal, Advances in Physics, 
will accept papers for publication in experimental and theoretical physics. 
The Philosophical Magazine publishes contributions describing new results, 
letters to the editor and book reviews. Advances in Physics publishes articles 
surveying the present state of knowledge in any branch of the science in which 
recent progress has been made. The editors welcome contributions from 
overseas as well as from the United Kingdom, and papers may be published 
in English, French and German. : 


acre ey 


The Magnetic Susceptibility and Electrical Resistivity of Some 
Transition Metal Silicides} 
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ABSTRACT 


Sintered samples have been prepared of a number of silicides formed by 
the Groups IVa, Va and VIa transition metals and silicon. Magnetic 
susceptibility measurements have been made at room temperature and 500°c 
and the electrical resistivity and temperature coefficient of resistivity 
determined for the range 20°c to 120°c. The results are discussed in relation 
to the type of bonding present in these compounds and it is shown that the 
transition metal silicides are essentially metallic in character. 


$1. [INTRODUCTION 


THE transition metals form with silicon a series of compounds which are 
in many ways analogous to the transition metal carbides. As for carbides, 
the most stable silicides are those containing metals from Groups [Va 
Va and Via of the periodic table. Many of these silicides have high 
melting points and their apparent refractory nature is evidenced by 
extreme room temperature brittleness. The limited and sometimes 
conflicting published data on the electrical properties of these compounds 
show, however, that they possess high electrical conductivity and, at 
least in the case of molybdenum disilicide, a positive temperature 
coefficient of resistivity (Glaser 1951). The present paper reports the 
results of some magnetic and electrical measurements which have been 
carried out on some of the more stable silicides in order to gain information 
on the type of bonding present in these materials. 


§ 2. EXPERIMENTAL METHOD 


It has been shown that specific homogeneous silicides can be prepared 
by mixing together silicon powder with the required quantity of metal 
powder and then heating a compact of the mixture to produce reaction 
(Robins and Jenkins 1955). This method of preparation has been 
employed in the present investigation using the same starting materials 
as in the earlier work. The mixed powders were pressed into bars and 
then heated in a static atmosphere of argon (99:7%) to a temperature 
at which reaction took place. The temperature required varied between 
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1300°c and 1550°o for the different mixtures. The reaction product 
was crushed to a powder which would pass through a 200 mesh screen 
and the powder was subjected to x-ray examination to confirm that the 
material was both single phase and homogeneous. 


2.1. Magnetic Susceptibility Measurements 


Susceptibility measurements were made on both powdered silicides 
and on sintered compacts by measuring the force on a sample held in an 
inhomogeneous magnetic field. The sintered compacts were prepared 
by compacting powder in a brass die at 10 tons/in? and then treating 
the compacts in the manner listed in table 1. It had been shown pre- 
viously that no change in the silicide would be expected to take place 
during the sintering process (Robin and Jenkins 1956). After sintering, 
the compacts were all ground to the same size (6.30 mm dia, and 1-25 mm 
thick) care being taken to prevent any possible contamination by 
ferromagnetic impurities. 

For the examination of powders, the samples were placed in a small 
flat-bottomed silica container (6 mm dia.) which was filled to a constant 
depth of 15mm with a known weight of powder. The container was 
supported by a silica rod between the flat pole pieces of an electromagnet, 
and the upper end of the silica rod was attached to a Sucksmith balance. 
For measurements on sintered samples the specimens were attached to the 
silica rod by a small cradle of thin platinum wire. During the experi- 


ments at room temperature the pole faces of the magnet were spaced - 


1 cm apart, but for measurements at elevated temperatures the separation 
was increased to 2-5cm to allow a non-inductively wound platinum 
furnace to be placed around the specimen. Measurements on each 
sample were carried out at three different current settings for the magnet 
and at least three determinations of the displacement were made at each 
setting. Before each measurement the electromagnet was moved in a 
vertical direction until the specimen was in that part of the field where 
H (dH /dx) was largest, this position giving the maximum deflection. The 
measured deflections were corrected for the effect of the holder and 
supports. 

Since the method of measurement of susceptibility employed was only 
comparative, a sample of high purity tantalum ground to the same 
dimensions as the sintered compacts was used as a standard sample. 


To confirm that the measuring technique was satisfactory, titanium and — 


zirconium compacts were prepared by compacting and sintering metal 
hydride powder in a vacuum (Robins 1956) and the measured suscepti- 
bility values were compared with the published values for these two 
metals. 

All the measurements on powders were only carried out at room 
temperature but, for many of the sintered samples, the variation in 
suceptibility over the temperatures range room temperature to 500°¢ 
was determined. The measurements were carried out in air and, to 
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confirm that significant contamination by oxygen or nitrogen at the 
elevated temperatures had not taken place, the room temperature 
susceptibility measurements were made both before and after the high 
temperature measurements. 


2.2. Hlectrical Measurements 


For the resistivity measurements bars of silicide, approximately 
11 cm long, | cm wide and 0-5 cm thick were prepared by pressing the 
silicide powder in a steel die at 40 tons/in® and then sintering the pressed 
bars using the treatments listed in table 1. The density of the sintered 
bars was measured and compared with pycnometrically determined 
values using the —200 mesh powder. The percentage of residual 
porosity could then be calculated. Densities of the silicides calculated 
from published x-ray data where these are available are also included 
in table 1. A considerable amount of residual porosity was present in 
some of the sintered bars and in order to assess the effect of this porosity 
on both the room temperature resistivity and the temperature coefficient 
of resistivity a number of bars of molybdenum disilicide and tungsten 
disilicide were prepared under different conditions to produce material 
of different densities. Only a limited amount of experimental work was 
carried out with hot pressed material due to the difficulties of preventing 
contamination by oxygen, nitrogen or carbon when using this technique 
with fine, reactive powders. 

The resistivity measurements were carried out using a Kelvin bridge 
having an accuracy of better than 0-5°,, and for both room temperature 
and elevated temperature measurements the bar and supporting frame- 
work were contained in a long tube furnace wound non-inductively. 
Resitance measurements were mace up to a temperature of 160°C and the 
temperature coefficient of resistivity calculated for the temperature 
range 20°C to 120°c. 


§ 3. MAGNETIC SUSCEPTIBILITY RESULTS 


Assuming the room temperature magnetic susceptibility of tantalum 
to be 0-84 10-® e.m.u./g (Kriessmann and Callen 1954), average values 
for the susceptibility of sintered samples of titanium and zirconium 
were found to be 3-24 10-® and 1-25 10-% e.m.u./g, respectively. A 
comparison of these results with the accepted values for these elements 
(3-2 10° °e.m.u./g for titanium and 1:3 10°%e.m.u./g for zirconium) 
showed that the experimental technique was satisfactory. The variation 
of the susceptibility of tantalum with temperature over the range room 
temperature to 500°C also agreed with the results of Kriessmann and 
Callen (1954), the paramagetism decreasing linearly with increasing 
temperature from 0-84 10~® to 0-78 x 10-6 e.m.u./g. 

The results of the susceptibility measurements on both powdered and 
sintered silicide samples are given in table 2 each value being the average 
of at least three readings. Paramagnetic values are given as positive 
and diamagnetic values as negative. In most cases results are given for 
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determinations at three different currents in the electromagnet but, 
particularly when carrying out measurements with powders having a small 
susceptibility, the deflections produced at the lowest current were too small 
for accurate measurement. In each case the value obtained at 12 amps 
was the most accurate, due to the larger deflection, and the accuracy of 
these values is considered to be better than --0-05 x 10-e.m.u./g. All the 
values have been corrected for the effect of the holder, this correction 
being small except where the sample has a small susceptibility. For solid 
specimens the deflection due to the platinum wire holder corresponded 
to approximately 2° of the deflection produced by the calibrating 
sample, and for the holder used with the powdered samples to 
approximately 5°,. 

No marked variation of susceptibility with field strength was found for 
any of the samples showing that ferromagnetic impurities were not 
affecting the results. It was known, however, that some of the powders 
used in the preparation of the silicides contained traces of iron and the 
absence of ferromagnetic impurity in the silicide was possibly due to the 
ferromagnetic element being present in the combined state. To confirm 
this a quantity of molybdenum disilicide powder was ground in a steel 
ball-mill to produce iron contamination, and subsequent magnetic 
measurements on this powder showed that significant contamination 
had taken place. Sintered compacts were then prepared from this 
contaminated powder and their susceptibility redetermined at different 
field strengths. The sintering treatment had caused all trace of the 
ferromagnetic impurity to disappear. 

Table 2 also contains susceptibility values for a number of the 
disilicides at 500°c, and for all the samples examined no marked deviation 
from linearity was found for the variation of susceptibility with temper- 
ature over the range room temperature to 500°c. High temperature 
susceptibility measurements were not made on the niobium disilicide 
and tantalum disilicide samples, due to the possibility of excessive 
contamination by oxygen caused by the reactive nature of these materials 
and the high porosity of the sintered compacts. 

The susceptibility values obtained at the maximum field strength 
for both room temperature and 500°c are recorded in table 2 as electro- 
magnetic units per gram molecule, and assuming that the value of 
— 3-65 x 10-® e.m.u./g atom for the susceptibility of silicon (Bates 1951) 
remains unchanged on combination with a transition metal, the magnetic 
susceptibility for a gram atom of the combined metal has been calculated. 
Any error introduced by this assumption is likely to be small since a 
change in the character of the silicon to Si*t(—1 x 10-e.m.u./g ion— 
Selwood 1956) would alter only slightly the value obtained for the 
susceptibility of the metal atoms. 

The disilicides of the first period transition elements were all para- 
magnetic, whereas the disilicides of the two later transition periods all 
exhibited diamagnetism. Increasing temperature increased the tendency 
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— towards diamagnetism, although the magnitude of the effect varied 

considerably for different compounds. In a number of cases the 
_ diamagnetic values obtained were numerically rather high. For the 

lower silicides of titanium and zirconium there was a general tendency 

for increasing silicon content to reduce the paramagnetism, although the 
_ titanium in the compound Ti,Si, had a slightly greater paramagnetism 
than unalloyed titanium. Conflicting data have been published regarding 
the compounds present in the zirconium-silicon system (Kieffer et al. 1954, 
Lundin et a/. 1953), but a compound corresponding to Zr,Si,; was prepared 
and its magnetic susceptibility was found to be +0-46 x 10-° e.m.u./g 
corresponding to +52 x 10°* e.m.u./g atom of metal. For the zirconium— 
silicon compounds, therefore, increasing silicon content did systematically 
reduce the paramagnetism of the metal. 

Only for chromium disilicide were there any previous data with 
which the present results could be compared. Foéx (1938) reported 
for this compound a weak, temperature independent diamagnetism 
(—0-094 x 10-® e.m.u./g), whereas the material examined in the present 
investigation possessed a weak paramagnetism. 


§ 4. ExecrricaL Resistivity RESULTS 


The effect of residual porosity on the electrical resistivity of sintered 
molybdenum disilicide and tungsten disilicide compacts is shown in 
figs. 1 and 2, and it is apparent that the porosity had a greater effect on the 
resistivity than could be accounted for by the presence of isolated pores 
assuming that isolated porosity only reduced the effective area of cross 
section of the sintered bars. Examination of the samples showed, 
however, that only isolated porosity was in fact present, the size of the 
pores decreasing with increasing sintering temperature. To confirm 
that the different sintering temperatures for the various bars were not 
producing a change in the basic material, additional bars of molybdenum 
disilicide were prepared varying both the sintering time and the sintering 
temperature to give material of the same density by different sintering 
treatments. The resistivity values again fell on the curve given in 
fig. 2. As additional confirmation, sintered material was crushed to 
powder, compacted and resintered and again the resistivity values 
corresponded to the curve in fig. 1. The reason for the large variation in 
resistivity with density for the molybdenum and tungsten disilicide 
samples was not established. It was found, however, that although this 
variation was large, no marked variation of the temperature coefficient 
of resistivity with density was observed. Temperature coefficient of 
resistivity values determined on bars containing an appreciable amount 
of porosity are, therefore, likely to be more reliable than the resistivity 
values ; these latter values are likely to be too high. 

Table 3 lists the room temperature resistivity and the temperature 
coefficient of resistivity over the temperature range 20°C to 120° for the 
compounds studied in this investigation, the resistivity values for 
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molybdenum and tungsten disilicides being taken from figs. | and 2. 
Published data on resistivity, where these are available, are also listed 
in table 3. The high conductivity of titanium disilicide as found by 
Glaser and Moskowitz (1953) was confirmed, its resistivity being less 
than that of the pure metal. The resistivity of the zirconium disilicide 
samples was a little lower than that found by Gallistl for hot pressed 


Fig. 1 


— 7 


| (a) Theoretical variation assuming | 
no interconnected porosity 


(b) Experimental values 


Resistivity (microhrns—cm) 


60 70 80 90 100 
Density (% of 6:10g/em?) 


The variation of room temperature resistivity with density for sintered 
molybdenum disilicide. 


material (and reported by Schwarzkopf and Kieffer 1953), but the 
resistivity was still more than double that of zirconium metal and con- 
siderably greater than the resistivity of titanium disilicide. The 
conductivity of the sintered vanadium disilicide was a little less than that 
reported for hot pressed material and this difference is possibly associated 
with a difference in residual porosity. A detailed study of the resistivity 
of niobium and tantalum disilicides was not carried out due to the large 
amount of porosity in the sintered sanyples, but some preliminary 
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measurements on tantalum disilicide gave a resistivity value considerably 
greater than that obtained by Gallistl as might be expected from the 
low density of the sintered material. The measured resistivity of 
_ molybdenum and tungsten disilicides agreed well with published data 
for these compounds (Glaser 1951, Schwarzkopf and Kieffer 1953). 


Fig. 2 


(a) Theoretical variation assuming 
no interconnected porosity 


(b) Experimental values 


50 


45 


Resistivity (microhms -—cm) 


60 7O 80 


Density (% of 9:62 g/m’) 


The variation of room temperature resistivity with density for sintered 
tungsten disilicide. 


Chromium disilicide was exceptional in having a very much higher 
resistivity than any of the other silicides examined. It is interesting to 
note that the most stable titanium-silicon compound Ti,Si, had a lower 
conductivity than either titanium or titanium disilicide. 

All the “compounds examined exhibited a positive temperature 
coefficient of resistivity and there was a general tendency for a large 
temperature coefficient to be associated with high room temperature 
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conductivity. Over the temperature range investigated the resistivity 
increased linearly with temperature, but the temperature coefficient for 
molybdenum disilicide of 6-38 x 10°* was considerably greater than that 
of 1:3 10-3 calculated from data given by Glaser (1951) for the temper- 
ature range 22°c to 65°c. 


Table 3. Electrical Properties 


Resistivity at 20°C 
(microhms em) Temperature 
Compound coefficient 
Present Published (« 10) 
investigation data 
TiSi, 16-7 123 4-63 
19 
ZrSi, 106-2 161 1-30 
VSi, 13-3 9-54) 3°51 
NbSsi, — 6:3) — 
TasSi, (38-0) 8-5) 3-32 
CrSi, 1420 > 250) 2-93 
MoSi, 21:8 21-5 6-38 
WSi, 38-2 33-44) 2°91 
Ti,Si, 350 — 0-86 
TiSi 39-3 — 4-13 
Zrsi 49-4 — 3°52 
| 


) Glaser 1951. 

(2) Glaser and Moskowitz 1953. 
(3) Kieffer e¢ al. 1953. 

‘) Schwarzkopf and Kieffer 1953. 


§ 5. Discussion 
In order to facilitate discussion of the results, the data obtained in this 
investigation are summarized in table 4 ; the room temperature resistivity 
values for the silicides of vanadium, niobium and tantalum are those 
given by Gallistl (and reported by Schwarzkopf and Kieffer 1953). 
Table 4 also lists for comparison the room temperature resistivity, 
temperature coefficient of resistivity, and magnetic susceptibility of the 


transition metals. The variation of magnetic susceptibility with temper-_ 


ature exhibited by the transition metals is in all cases small, each element 
in a particular group showing the same type of temperature dependence ; 
the sign of the temperature coefficient alternates from one group to the 
next, the paramagnetic susceptibility increasing with temperature for the 
Group [Va and Vla elements and decreasing with temperature for the 
Group Va elements (Kriessman and Callen 1954). 

The following features of the experimental results may be put forward 
as evidence for the metallic character of the transition metal silicides. 
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1. The electrical conductivity of the silicides is high and in a number 
of cases is, in fact, greater than the conductivity of the uncombined metal. 

2. The temperature coefficient of resistivity is in all cases positive 
and of the same order as for the transition metals. 

3. The silicides of the metals of the first transition series exhibit 
paramagnetism, the magnitude of which may be associated with “free ° 
electrons and there is no reason to expect a marked difference in the 
nature of the silicides of the first and later periods. 

4. The variation of the susceptibility with temperature is small. [f 
the Weiss equation 


m= CT 8) i» b> 


is assumed to hold then the results for titanium disilicide and vanadium 
disilicide, yield values for the Curie temperature (@) of —3050°K and 
—720°K respectively, and Bohr magneton values, calculated from the 
Curie constant (C), of 1-9 and 1-2, respectively. The large negative 
values of 6 would suggest that the susceptibility cannot be interpreted in 
terms of the classical theory of paramagnetism and that any ionic contri- 
bution must be small. 

5. The magnetic susceptibility of the various compounds in the 
titanium-silicon and zirconium-silicon systems changes in an almost 
symmetrical manner with increasing silicon content. 

It is apparent, therefore, that these compounds are essentially metallic 
and they should be considered either in terms of band theory or the more 
empirical resonating covalent bond approach. 

When considering the experimental results in more detail it must be 
borne in mind that the disilicides do not all possess the same crystal 
structure (Robins and Jenkins 1955) and in each case the structure is 
unrelated to that of the original metal. An analysis of the various 
interatomic distances has shown that the metal-metal distances in the 
disilicides are so large that any significant amount of bonding between 
the metal atoms seems improbable and, therefore, the metallic behaviour 
of these materials must be associated with the character of the strong 
metal-silicon bonds (Robins and Jenkins 1955). 

There is a general correspondence between the pattern of the magnetic 
susceptibility results for the disilicides and that of the transition metals, 
although all the transition metals exhibit paramagnetism, whereas many 
of the silicides are diamagnetic. For both the metals and the silicides, 
the tendency towards diamagnetism increases with increasing atomic 
number in a given group, although only the silicides in the first transition 
series are paramagnetic. The transition metals exhibit alternately 
relatively high and low values of room temperature susceptibility for the 
first five elements in a given transition series (Stoner 1954): the Group 
Va elements have high values and the elements of Group [Va and Via 
relatively low values. This is also true for the disilicides. Vanadium 
disilicide exhibits the highest paramagnetism of the various compounds 


Resistivity of some Transition Metal Silicides 325 


examined and the disilicides having numerically the highest diamagnetism 
contain elements from Groups [Va and VIa. 

The large paramagnetism of the transition metals is due to the high 
density of states at the top of the Fermi distribution, resulting from the 
interaction of incomplete d shells. Stoner (1954) has shown, however, 
that estimated values of the magnetic susceptibility of these elements 
calculated from the electronic contribution to the specific heat, assuming 
any additional interaction to be negligible, are considerably less than the 
observed values. There must be present in the transition metals, there- 
fore, exhange interaction contributing to the observed paramagnetism. 
Although no direct comparison can be made between the magnetic 
properties of the silicides and those of the transition metals, it would 
seem probable that the increased metal—metal separation in the silicides 
would reduce both the density of states at the top of the Fermi distri- 
bution and also the magnitude of any exchange forces between the 
incomplete d shells. However, the value of the susceptibility of titanium 
calculated from the electronic specific heat is given by Stoner (1954) as 
48x 10-%e.m.u./g atom which is lower than the susceptibility after 
combination with silicon. In titanium disilicide, therefore, in addition 
to a high density of states there must also be a significant contribution to 
the paramagnetism from exchange forces. It would seem likely that, 
for vanadium disilicide also, the high paramagnetism can only be 
explained by assuming both a high density of electron states at the top 
of the Fermi distribution and a significant contribution from exchange 
forces. Since, as has been shown above, the metal-metal distances are 
such that direct interaction between incomplete d shells seems improbable, 
the high paramagnetism of titanium and vanadium disilicides lends 
support to Zener’s hypothesis (1951 a,b, 1952) that the d shells are 
indirectly coupled through the conduction electrons. 

The high paramagnetism of titanium disilicide and vanadium disilicide 
is associated with a high electrical conductivity which in both cases is 
considerably greater than the conductivity of the unalloyed metal. Since 
in the transition metals the high electrical resistivity is caused by the 
large scattering effect of the high density of states in the incomplete d 
bands, the low resistivity of these silicides would suggest that the scattering 
effect is reduced due to a decrease in the density of states. On Pauling’s 
hypothesis this increase in conductivity may be associated with the role 
of the silicon atoms as centres of pivoting resonance thereby reducing the 
contribution from the scattering effect of the metal atoms. 

The marked difference between both the electrical and magnetic 
properties of titanium disilicide and zirconium disilicide is most striking, 
particularly since the crystal structures are very similar. The high dia- 
magnetism of zirconium disilicide would suggest that the zone structure 
is such that the effective electron mass of the electrons in the highest zone 
is small. This interpretation of slight overlap is supported by the high 
electrical resistivity of this compound which is an order of magnitude 
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greater than that of titanium disilicide. It is possibly of significance that 
the interatomic distances in zirconium disilicide are greater than in 
titanium disilicide. 

The high electrical conductivity of niobium disilicide and tantalum 
disilicide would suggest the presence of ‘free’ electrons which would 


Ss 


produce a weak paramagnetic effect and the observed weak diamagnetism — 


is, no doubt, due to the diamagnetism of the completed shells and paired 
electrons being greater than the Pauli paramagnetism. The absence of a 
large paramagnetic effect as is present in vanadium disilicide shows that 
any exchange interaction must be small possibly due to the interatomic 
distances in the silicides increasing with increasing atomic weight of the 
metal. 

The results obtained for the disilicides of the Group Vla metals present 
a number of interesting features. Chromium disilicide, unlike the 
disilicides of molybdenum and tungsten, has the same crystal structure 
as the Group Va disilicides. It exhibits a weak, almost temperature in- 
dependent, paramagnetism which is coupled with a very high electrical 
resistivity. Molybdenum and tungsten disilicides, however, both have 
comparatively low resistivity values, and are diamagnetic. The numerical 
value of the diamagnetism of molybdenum disilicide increases considerably 
with temperature as also does the electrical resistivity. A theoretical 
analysis of the band structure of molybdenum and tungsten disilicides has 
been carried out by Schenk and Dehlinger (1956) who predict from their 
results that these compounds should possess diamagnetism. As for 
niobium disilicide and tantalum disilicide, the Pauli paramagnetism must 
be more than balanced by the underlying diamagnetism. 


§ 6. CONCLUSIONS 


The high electrical conductivity of the silicides studied in this investi- 
gation and their positive temperature coefficient of resistivity showed that 
these compounds are essentially metallic. 

The disilicides of the elements in the first transition series are para- 
magnetic and for titanium disilicide and vanadium disilicide, the magni- 
tude of the paramagnetism is comparable with that of the transition 
metals. This would suggest that for those compounds, in addition to a 
high density of states at the top of the Fermi distribution, there must also 
be a significant contribution to the paramagnetism from exchange forces 
in spite of the large metal-metal separation. The disilicides of the 
elements in the second and third transition series are all diamagnetic, 
zirconium having a particularly large diamagnetism coupled with a 
relatively high electrical resistivity. The variation of magnetic sus- 
ceptibility with temperature over the range room temperature to 500°C 
was small, for the compounds studied. 

It is apparent from the results that the study of this type of compound 
may be of importance as a method of obtaining additional information on 
the nature of the transition metals, 
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ABSTRACT t 

The process 7+ > 7°+et+t-+yv has not yet been observed. Estimates of its 
branching ratio relative 10 7+ — ++ can be derived (1): from a phase space 
volume argument in conjunction with the available experimental information 
on wt and (assumed analogous) K+ processes, and (2): from a calculation 
of the rate of the beta decay transition between the ‘z+’ and ‘7°’ states 
of a bound nucleon—antinucleon structure. 

A possible experimental scheme for detection of a+ —>7°--et-+-y in the 
presence of the overwhelmingly more frequent 7+ + ++» is outlined. 


§1. INTRODUCTION 


THE process 7+ > 7°+e++v has not up to now been observed. This 
process is interesting in the context of known beta decay modes of 7 and 
kk particles. Empirical information on corresponding 7 and K* beta 
decay modes is exhibited in the table. The experimental evidence on 
a,b, d,and a’ is consistent with the interpretation of the neutral unobserved 
particles as a neutrino. 
We first note that processes 6 and d are about equally probable; hence 
each of the two reactions 
K= > 7° + pF+p oily gh, em 


K+ n+et+p Lye ts 


may have the same type of coupling and the same coupling constant. 
(‘Lhe phase space volume factor in the transition probability for the muon 
mode 6 is only slightly smaller than for the electron mode d.) 

Modes ¢ and c’ have not so far been observed. In particular the 
branching ratio of 7* decay by c’ is now known to be less than 10~*8. 

The direct analogue of b is b’; however b’ is energetically forbidden and 
can occur, if at all, only as a virtual process. 


+ Communicated by the Authors. This research was supported in part by 
the Office of Naval Research of the Navy Department and the Air Force 
Office of Scientific Research of the Air Research and Development Command. 

{ A preliminary version of this note has been given by Feenberg and Primakoff 
(1957). See also the discussion by Weinberg (1957). ‘ 

§ The most recent experimental results on c’ are due to Anderson (1957). 
The branching ratio of 7+ decay for the radiative form of c’: a+ > Chaya 
is (3+7)x10~® (Cassels 1957), 
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These remarks suggest the possibility of a close analogy between each 
of the 7 decay processes a’, ...d’ and the corresponding K processes 
a,...d. We use this analogy to obtain an estimate of the rate of the 
as yet unobserved process d’: 7+ 7°+et+v. The argument leading 
to our estimate is fairly arbitrary so that we consider that the result has 


Information on Corresponding A+ and 7+ Decay Processes} 


Beg 
P (ut or et)t | 
MH 66 Partial 
Process @Q/m .c? or decay 
Pmax Ae el ae rate 
Pats (ut or e¥)t 
(a) K+ + wt+p 759 462 5 X 107/see 
(b) K+ + 7®+pt++yv 495 422 4 108/see 
(c) Kt + et+typ 965 483 not observed 
(d) Kt + w®+et+y 701 448 4x 108/see 
(a')  wt+>ptty 66 57 4x 107/see 
, “4. nfs eee 2" Re energetically 
al iret ae ee forbidden 
(c’') at—setty 272 136 not observed 
(’)  at+n®*+et+y 8 9 q 


+ Shapiro (1956). + Computed. 


value only as a guide in the search for reaction d’. Clearly, observation 
of this reaction would be of considerable fundamental interest, but even 
experimental evidence that its branching ratio is appreciably smaller 
than our estimate would be helpful in the development of a coherent theory. 


§2. ESTIMATE OF BRANCHING Ratio 


Let now M_, denote the effective transition matrix element for process 
a(x=a,b,...,d’), p, the associated number of final states per unit final 
energy range, i.e. the associated ‘ phase space volume’, and #, the corres- 
ponding partial decay rate. Then, with p,, p,, p, representing the 
neutrino, muon, and electron momenta, we have: 


2 
R,= > |M.|ex ie eRe a he ee hs (1) 


0b 
_ 2 spat 
Px = (const) p, / ( an 


r PDmax r 


a | OE s 
= t dp p,2 | dé sin 6,472 (S")} =P, OC Pes 
Pz = (cons ) | 5 PP} ix ty SIN U), \P» j ap, Li P; 
x=b,d,d’. (3) 


Jie=a case Ret een) 


2 


P.M. 2A 
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Thus, writing m,+,m,+,M,0, m, for the kayon, charged pion, neutral 
pion, and muon masses, and with use of the table, we have the processes 
a’ and a: 


[Ma |? pv s, eee m,+* mM, ; )/m+ [1 + (m,/m,+)?] ' (4) 


te rf [wr Pa 2 |M, | 9 M K Re a mM,, \/mg+ [1 + +(m m_,|m«+)?] 
4x 10’ 
~ Bx 10? 
Consequently, 


ee, as 


a 


The corresponding calculation for the processes d’ and d yields: 


Ke a | Ma Pew ~ | Mw |? ee TS ee a = | Mal? 1.9 x LO) 
Ra [Mal? pa |Ma|? (mg+)?/270 |u|? 


Ry 


Tx 108 


We now assume that the ratio of corresponding 7 and K effective transi- 
tion matrix elements is independent of the decay channel, i.e. using also. 


eqn. (5): 


M V MY My My ’ 
“| ~ ~ |-—!|~ aaMiy A 6 See 
M,| \m\ \™, : 7) 
Eqns. (6), (7) and the table then give: 
Es 2 O-2180C st gl tige oe kOe. ie = ee 


so that, according to eqn. (8), DIMES out of every 10° charged. 
pions are expected to decay via m+ —> 7° +e++v. 

An alternative estimate of R, is obtained by interpreting process d’ 
as a beta decay transition analogous to the beta decay of a bound system 
of nucleons and characterized by the same Fermi and Gamow-—Teller 
types of interaction with the same coupling constants. Introducing the 
corresponding effective ‘nuclear’ matrix element N, (Nz proportional 
to M,) and the f function of beta decay theory, and interpreting the 
transition as superallowed within a 7'=1 isobaric spin multiplet, we 
obtain : 

|N @ [27(9m 2) 2 ~ 5000 ee 
a’ 
with 
[Ma |2= |(P=1, 1)=0|1,-i7,|7=1, T,=1) P= 


Also f(9m,c?)~1700. Equation (9) then yields 
Ra = 0-5/sec ; RaylRy& 10-8 co 5 5 a é, (10) 


quite close to the estimate of eqn. (8). 
We may now attempt to place the above calculation on a possibly more 
secure basis by using a bound nucleon—antinucleon (Fermi-Yang) model 
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for the pion. In the two particle approximation of this model the r+, 7° 
wavefunctions arey : 


WS [p(1)n*""(2) Jr, Pelv roti Xa Sieg) ee es CEE) 
m= 2-N2([ (1) (2)] — [(1)n*™(2)]) naa, 7, <0 7—0-(1, Xe 81, 89). 
In this context the transition occurs via one or the other of the two processes 


(P)bouna a (2) round ae et ald 2 : < - G ( 12), 
‘Guaa eee = A Saige Pee aa en + v. 


These processes may contribute in-phase or 180°-out-of-phase to the: 
a*—>7® transition amplitude depending on the type of coupling between 
the nucleon—antinucleon and lepton fieldst. In the in-phase case— 
provided that the coupling is such that the 7+->7° transition is not 
otherwise * forbidden’ {—the estimate of eqns. (9) and (10) for the branching 
ratio R,/R, is valid, but in the out-of-phase case the computed R,/R,,- 
will be much smaller than 10-8. 


§3. PossIBLE EXPERIMENTAL SCHEME 


We shall now discuss briefly a possible arrangement for detecting the 
process 7*—>7°+e*+v in the presence of an overwhelming background 
from the process 7*->y++v. <A beam of positive pions is stopped in an 
absorber. The process 7+->7°+e++v then gives rise to two pairs of 
photons from the reactions: 


0 
7 > Y132m ec? + Yis2mect 
+ —-. 
Sip lor Ymec? A VY mec? 


with the m,c? photons from the e*, e~ annihilation delayed relative to the 
132m,c? photons from the 7° decay by less than 10-* sec. On the other 
hand, the process 7++pt+v is followed by 


pot—>er+ v+v; et+ CVn ect T Vinee? 


with any accompanying high energy (<104m,c?) photons arising only 
from external or internal bremsstrahlung of the decay e+ from the pt. 


+ In eqn. (11), the %)~0-(xy, X2, $1, Sg) are the angular momentum, parity 
eigenfunction factors in the 7+, 7° wavefunctions, x,, X», $1, $2 being the space 
and spin coordinates of the nucleon and antinucleon. For the forms of the 
nucleon—antinucleon T?, 7', eigenfunction factors in eqn. (11), see for example, 
Malenka and Primakoff (1957). 

t Thus, with the now current beta decay interaction Hamiltonian (which 
includes any existing non-invariant terms for charge conjugation, space reflection 
and time inversion in the lepton factor) the two processes (12) contribute in-phase 
for polar vector or tensor coupling and contribute 180°-out-of-phase for scalar, 
axial vector or pseudoscalar coupling. However, since the 7+ > 7° transforma- 
tion involves a J —0- + J=O- transition, the tensor coupling contribution for 
each of the two processes (12) and so for their (in-phase) combination is 
‘second forbidden ’ and hence is much smaller (for given magnitudes of the 
coupling constants) than the ‘ allowed’ polar vector coupling contribution. 


2A2 
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With reference to fig. 1, the approximately oppositely directed 132m ,c? | 
photons can be detected in coincidence. Energy discrimination in the 
counting circuit can be used to decrease the number of background 
coincidences arising from, e.g. the above mentioned high energy brems-— 
strahlung photons accompanying the »+—e* decay. If the production — 


‘ 
‘ 
4 


Fig. 1 
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(32meC)¥ 
ABSORBER 


COINCIDENCE AND 
DISCRIMINATOR CIRCUIT 


Block diagram of apparatus for detection of neutral pions produced in absorber. 


of neutral pions in the absorber is demonstrated by this arrangement, an 
energy discriminating delayed fourfold coincidence circuit for the detection 
of the two 132m,c? photons and the two m,c* photons—see fig. 2—can be 
used to determine whether the neutral pions are indeed produced by the 
process 7+—+>7°+e++v or whether they arise from a charge exchange 
reaction in the absorber : 


at++(Z,A)>7°+(Z+1, A)f. 
Another possible source of background will arise from the conceivably 


existing radiative 7+-—>e* decay followed by (external) bremsstrahlung 
of the e+ and its ultimate annihilation : 


art > er+tut+y; et—>(et)'+y'; (6°) Fes ¥,,, at Vneee 


+ The delay involves the time interval dt between the instant of passage 
of the 7+ through a monitor counter in front of the absorber and the 
instant of detection of the (energy discriminated) fourfold coincidence. The 
number of delayed fourfold coincidences arising from the process sought : 
at > 7+ et+-+-y, ete., will vary with At as exp (— At/r) (r=7* lifetime —2-5+ 108 
sec). On the other hand, the number of delayed fourfold coincidences associated 
with neutral pions which are not decay products of incident positive pions, 
will exhibit a dependence on At quite different from exp (—At/r). Thus, the 
delayed fourfold coincidence circuit should discriminate effectively against 
the 132m,c? photons from any promptly produced charge exchange neutral 
pions—there will also be effective discrimination, on a time basis as well as on 
the previously noted energy basis, against the relatively late appearing high 
energy bremsstrahlung photons associated with 

at —>pt+y; pt—>ettytyty; et (e*)'-y’. 
Finally, the number of prompt charge exchange neutral pions can be kept to @ 
minimum by use of incident positive pions of low energy falling on a high 
Z absorber while the probability of external bremsstrahlung by the decay e+ 
from the «+ can be decreased by making the absorber as thin as feasible. 
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In general however, the two photons y, y’ will not be emitted mutually 
antiparallel and will not each have an energy close to 132m,c2. On the 
other hand, one may search for a+->et++ p+ VY; Cte > Ym ett Vineet 
with an appropriate modification of the apparatus of fig. 2. . , 


Fig. 2 
PION MONITOR COUNTER 


COUNTER COUNTER 
‘BEAM 


ou 
: Yoe 
Ns) 
TWOFOLD COINCIDENCE ~ — TWOFOLD COINCIDENCE 


AND DISCRIMINATOR AND DISCRIMINATOR 
CIRCUIT ay (m CIRCUIT 


g = 
COUNTER oo ee reer “S09 COUNTER 


DELAYED FOURFOLD COINCIDENCE CIRCUIT 


| DELAY LINE 


Block diagram of apparatus for detection of delayed four-fold coincidences. 
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ABSTRACT 


A new method is described for the decoration of dislocations in crystals 
of silver chloride. The crystals are sensitized by heating them, together 
with gold wire and cupric chloride, in an atmosphere of chlorine in sealed 
tubes at a temperature of 350°c. They are then strained, exposed to light 
with a wavelength which is not strongly absorbed by silver chloride, and 
heated for a few minutes at a temperature between 140 and 160°c. Nuclei 
which are formed along the dislocation lines during expesure increase to 
visible size during the subsequent thermal treatment. 

The formation of systems of prismatic dislocations around precipitates 
in silver chloride crystals is described and it is suggested that the observations 
may be relevant to the interpretation of precipitation hardening in certain 
alloys. 


§ 1. IyrRopUCTION 


DisLocations in crystals of silver chloride are made visible by exposing 
them to light, which causes silver to separate along the dislocation lines 
(Hedges and Mitchell 1953, Clark and Mitchell 1956, Mitchell 1957 c), 
and by etching the surfaces with a 3N solution of sodium thiosulphate 
(Jones and Mitchell 1957). Both the dislocations which remain after 
annealing and those which are introduced by plastic deformation are 
rendered visible by these two methods. The decoration of the disloca- 
tions in crystals of silver chloride and of silver bromide with photolytie 
silver does not depend on the introduction of foreign substances or upon 
thermal treatment of elevated temperatures. The configurations of 
dislocations which are observed in annealed and plastically deformed silver 
halide crystals are therefore likely to have been modified to a smaller 
extent by the decoration process than those in alkali halide crystals. 
The decoration of the dislocations is practically continuous and discrete 
particles of photolytic silver cannot usually be resolved. 

The only real disadvantage of this method is that it has never been 
possible to decorate dislocations at a depth greater than 30, below the 
surface, even when the crystals have been sensitized by the deposition of 
thin films of silver or of silver sulphide upon their surfaces before exposure 
to light. The aim of the present work has therefore been to discover 
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a method for making dislocations visible throughout the thickness (usually 
200 to 400) of the thin silver chloride sheets which we use for much of 
our work on plastic deformation. Our object has been achieved with 
S eeaprent of a new method for decorating dislocations with platelets 
of gold. 


§ 2. EXPERIMENTAL MerHops AND OBSERVATIONS 


During their work on the photographic sensitivity of mixed crystals of 
aurous and silver halides, Clark and Mitchell (1956) observed that thin 
platelets of gold, with surfaces parallel to {111} planes, separated 
spontaneously when the crystals were stored at room temperature after 
exposure to visible light. We hoped that it might be possible to find 
conditions in which the gold would separate upon nuclei formed along 
dislocation lines. 

The first stage in the present work was therefore the preparation of 
crystals of silver chloride with a range of concentrations of aurous chloride 
extending from 10+ to 10-3mol%. 10-'mol®% of anhydrous auric 
chloride, which decomposes to give aurous chloride, was added to the 
molten silver chloride and globules and thin sheet crystals were made as 
described by Clark and Mitchell. Hydrogen chloride was used as a 
protective gas throughout the operations. After separation from the glass 
plates, the crystals were cut into sections and annealed in a sealed tube 
in an atmosphere of chlorine at a pressure of 200 to 300mm Hg at 350°C 
for: 6 to 12 hours. 

The crystals were then exposed to the light from a 250 watt mercury 
vapour lamp which was filtered through a thin sheet of silver chloride for 
a number of the experiments to remove strongly absorbed wavelengths. 
After exposure, the crystals were heated in air on thin microscope cover 
slips resting on the surface of a heated steel bar. 

When annealed crystals containing 10-'mol% of aurous chloride were 
exposed and either kept at room temperature for several weeks or heated 
for a few minutes to temperatures between 140 and 160°c, they darkened in 
the exposed areas due to the separation of platelets of gold throughout 
their thicknesses; no precipitation occurred in unexposed areas. Indi- 
vidual dislocations could not be recognized in these crystals but the density 
of gold platelets was appreciably greater on grain boundaries and sub- 
boundaries than elsewhere. All the earlier observations of Clark and 
Mitchell were confirmed. 

Globules with progressively lower concentrations of aurous chloride were 
next prepared by diluting the 0-1 mol% material with silver chloride. Thin 
sheet crystals were made, cut into sections, annealed in chlorine and exposed 
and heated as just described. The density of gold platelets decreased with 
the concentration of aurous chloride and clusters in which platelets formed 
on all the {111} planes passing through a nucleus did not appear. The 
individual platelets of gold in a erystal containing 0-05mol% aurous 
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chloride are shown in fig. 1, Pl. 9. This crystal was exposed at a distance. 
of 25cm from the mercury vapour lamp for 6 minutes without the produc- 
tion of any microscopically visible internal change. It darkened rapidly 
in the exposed area due to the formation of platelets of gold when it was. 
heated for 1 minute at 140°c. These randomly distributed platelets, 
which have surfaces parallel to the {111} planes, are not formed when 
crystals with lower concentrations of aurous chloride are treated in the- 
same way. ; 

Crystals with 0-005mol% aurous chloride were then prepared and 
annealed in chlorine at 350°c for 8 hours. This produces recrystallized 
areas around the cut edges in which very few dislocations remain and 
polygonizes the central areas of the crystal sections. Gold platelets did 
not separate near the surfaces when these crystals were heated after 
exposure. Precipitates were, however, formed within the crystals at 
distances greater than 50 ,. from the surfaces and there was some decoration 
of dislocations. No precipitation was observed in crystals containing 
0-001 mol®%, aurous chloride. 

It seemed likely that some of the aurous chloride had diffused out when 
these crystals were annealed in chlorine at 350°c. We therefore heated a 
further batch of the crystals in a sealed tube with chlorine, together with a 
piece of gold wire, for 8 hours at the same pressure and temperature. The 
final concentration of aurous chloride after this treatment is not known but 
new and interesting phenomena were observed when the crystals were 
examined with the microscope after having been exposed for 1 to 5 minutes 
and heated for 5 to 15 minutes at temperatures between 140 and 160°c. 
After they had been annealed in chlorine, together with the gold wire, the 
crystals which had an initial concentration of 0-005 mol% aurous chloride 
contained microscopically visible internal particles, with a diameter of 
less than ly, distributed throughout the thickness. When they were 
exposed for 5 minutes and heated for 15 minutes at 160°C, arrays of 
regularly spaced smaller precipitates appeared along the twelve (110) 
directions passing through the centres of the 1, particles. This new 
precipitation phenomenon which is illustrated in fig. 2, Pl. 9 will be 
discussed in § 3. 

The dislocations of internal sub-boundaries in crystals which had an 
initial concentration of 0-001 mol% aurous chloride were decorated when 
they were exposed for 1 minute and heated for 6 minutes at 160°c. A 
decorated hexagonal network in the interior of a crystal is shown in fig. 3, 
Pl. 10; the high density of randomly distributed particles should be: 
noted. The smaller secondary precipitates which are shown in fig. 2, 
Pl. 9 did not appear in these specimens. All the thin sheet crystals. 
which were prepared from globules of silver chloride containing aurous 
chloride had a high density of particles after annealing and, although it 
proved possible to decorate the dislocations of the tilt and twist sub- 
boundaries, it was evident that such crystals would not be suitable for: 
work on the behaviour of dislocations during plastic deformation. 
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Thin sheet crystals of very pure silver chloride were next prepared and 
annealed at 350°c for 8 hours in chlorine at a pressure of 200 to 300 mm Hg. 
They were then further annealed at 350°c for 2 hours in chlorine at a 
pressure of 30mm Hg together with a piece of thin gold wire. This is 
essentially the technique which is described in the footnote to the paper 
by Barber et al. (1957). Gold does not separate spontaneously along 
dislocation lines when the silver chloride crystals are cooled down to room 
temperature and their properties differ in this respect from those of sodium 
chloride crystals. The dislocations remaining after the two annealing 
treatments were clearly decorated when the crystals were exposed for 
5 minutes and then heated for 10 minutes at 140°c. During the heat 
treatment microscopically visible particles of gold never separated in 
unexposed areas. The quality of the decoration is shown in fig. 4, Pl. 10. 
The discrete nature of the gold particles is illustrated by the dislocations 
in the plane of the photomicrograph near the upper left hand corner of the 
field. All the phenomena which have already been described for crystals. 
of pure silver chloride and silver bromide in which the dislocations have 
been decorated with photolytic silver have now been observed in crystals 
of silver chloride in which the dislocations have been decorated with gold, 
with the advantage that the decoration extends throughout the crystals 
instead of being confined to the surface region. 

In an attempt to improve the decoration, the pressure of chlorine in the 
tube during the second anneal was next reduced from 30 first to 15 and 
then to 5mm Hg at 350°c. Gold did not separate in any form when the 
crystals from these batches were heated to 160°C after exposure. Luckey 
(1953, 1955) had previously observed that silver bromide was sensitized for 
the liberation of bromine during exposure to light at temperatures around 
100°c by the incorporation of cupric bromide. We therefore decided to. 
investigate the effect of the introduction of cupric chloride on the separa- 
tion of gold during the thermal treatment of exposed crystals of aurous. 
and silver chlorides. By evaporating the required volume of a dilute 
solution to dryness, we deposited 5x 10-°gm of cupric chloride on the 
internal surface of the glass tube which was to be used for the second 
anneal. Previously annealed crystals were then heated for 4 hours at 
400° in chlorine at a pressure of 15mm Hg, together with a piece of gold 
wire. They were exposed for 5 minutes and heated for 5 minutes at 
160°c. The dislocations were clearly decorated throughout the thickness 
of the crystals by the separation of discrete particles of gold during the 
thermal treatment after exposure and there was a negligible background 
density of randomly distributed gold particles. The most interesting 
observations were made in the recrystallized areas around the edges, 
which have a very low initial density of dislocations. Dislocations 
introduced into these areas by plastic deformation after the diffusion of 
aurous and cupric chlorides into them are shown in figs. 5 and 6, PI. 11. 
The decorating particles are larger and more widely separated than the 
particles of photolytic silver which form in dislocations near the surface 
during exposure, but optimum conditions may not yet have been established. 
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An example of the detail which can be observed is given in fig. 7, Pl. 10. 
In this specimen, the boundary separating the recrystallized area from the 
remainder of the crystal moved discontinuously during the first anneal in 
chlorine at 400°c and left behind a number of microscopic cavities on the 
surfaces corresponding with the arrest positions ; some cupric chloride may 
have separated in these during the second anneal. The cusps formed when 
dislocation lines gliding in the recrystallized area were held back by these 
imperfections are clearly shown in the photomicrograph. 

We have established that there is continuity between the sub-surface 
decoration of dislocations with photolytic silver and the internal decoration 
of dislocations with gold by exposing crystals to the full radiation from 
the mercury vapour lamp and then heating them to 150°c for 5 to 10 
minutes. The dislocations of many tilt boundaries could then be followed 
downwards from the surface where the decoration was produced mainly 
by photolytic silver to the interior where the decoration was due to the 
separation of discrete particles of gold. 


§ 3. DIscussION 


In this section, we shall discuss the mechanism of the decoration of 
dislocations in crystals of silver chloride with gold and the significance 
of the precipitation phenomenon which is illustrated in fig. 2, Pl. 9. 
After it had been established that dislocations could be decorated as shown 
in fig. 3, Pl. 10, our work was concentrated on eliminating the randomly 
distributed particles which were observed both after the diffusion of the 
aurous chloride into the crystals and after the subsequent exposure and 
heating. The density of these particles was reduced (1) when the speci- 
mens were recrystallized by annealing them for a long period in chlorine 
before the aurous chloride was diffused into them ; (2) when the concentra- 
tion of aurous chloride was decreased to the smallest value at which gold 
would precipitate and (3) when cupric chloride was diffused into the 
specimens at the same time as the aurous chloride. 

The Schottky defects which are present in significant concentrations in 
crystals of silver chloride at temperatures near the melting point aggregate 
when the crystals are cooled to lower temperatures, with the formation of a 
random distribution of small dislocation loops and internal cavities. 
During the processes of decoration, gold nuclei are formed at these 
imperfections and grow into visible particles. Our observations show 
that the dislocation loops and cavities are usually eliminated when the 
areas around the cut edges of the thin crystals recrystallize during the 
course of a long anneal in chlorine at 350°C. Nuclei cannot then be formed 
in regions of the crystals free from dislocations unless vacant chloride ion 
lattice sites are available which can condense together with vacant silver 
ion lattice sites to produce internal cavities in which gold atoms can 
aggregate. At temperatures below 200°c, the concentration of vacant 
chloride ion lattice sites in silver chloride crystals is negligible; nuclei 
can therefore only be formed in dislocations and other crystal imperfections 
where space for the separation of gold atoms is already available. 
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Nuclei consisting of groups of silver and gold atoms are formed on the 
surfaces of internal cavities and along dislocation lines when the crystals 
are exposed to light with a wavelength beyond the absorption edge of silver 
chloride. The nuclei adsorb gold or silver ions and become positively 
charged when they exceed a certain critical size, the compensating negative 
charge being provided by vacant silver ion lattice sites (see Mitchell 1957 b 
for a discussion of these properties). When the crystals are heated to a 
temperature between 140 and 160°6, electrons are transferred from the full 
band (with the production of positive holes) to the low lying acceptor levels 
in the forbidden band associated with the positively charged nuclei. Their 
positive charge is restored by the adsorption of further cations and further 
vacant silver ion lattice sites are created. The positive holes and vaeant 
silver ion lattice sites diffuse to the surface of the crystals where chlorine 
escapes. Copper chloride has a number of functions during sensitization, 
exposure and heat treatment. One of them appears to be that of increasing 
the efficiency with which chlorine molecules are liberated from the surfaces 
of the crystals during heat treatment at temperatures between 140 and 
160°c. The cupric ions introduce acceptor levels in the forbidden band 
near the top of the full band to which electrons may be transferred by 
thermal activation with the creation of positive holes. A fraction of these 
positive holes, which are repelled by the positively charged gold nuclei, 
are trapped at the surface in the form of Cl, -molecular ions. These 
then trap positive holes, liberated when electrons are transferred from the 
full band to the positively charged gold nuclei, and form chlorine molecules 
which escape (see Mitchell 1957a,b for a detailed discussion of this 
process). 

We turn now to the discussion of the significance of the fine precipitates 
which are shown in fig. 2, Pl. 9. It may be recalled that the central 
particles of the systems of precipitates were visible in the microscope when 
the crystals were removed from the tubes in which they had been heated 
with the gold wire in chlorine. The systems of regularly spaced fine 
precipitates extending along the (110) directions passing through the 
central particles first became visible after the crystals had been exposed 
to light at room temperature and heated to 160°c for 6 minutes. 

Although the phenomena are on a smaller scale, there is a close analogy 
between these systems of fine precipitates and the systems of positive 
prismatic dislocations which form during cooling around glass spheres 
imbedded in silver chloride crystals and which have recently been described 
by Jones and Mitchell (1958). Asa working hypothesis, which is consistent 
with all the observations which have been made, we propose that the 
central gold particle grows initially at the annealing temperature by a 
mechanism which depends on the diffusion of Schottky defects to its 
interface with the silver halide. Precipitation continues during cooling 
in the furnace but a stage is reached at which the rate of arrival of pairs 
of lattice vacancies at the interface is insufficient to make space available 
for the gold atoms separating there. The gold nuclei then increase in 
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size until the decrease in free energy accompanying the separation of gold 
atoms is balanced by the increase in free energy due to elastic deformation 
of the surrounding silver halide crystal. A radially symmetric compressive 
stress field is thus established around the gold particles which is of the same 
nature as the stress field around the glass spheres in the experiments of 
Jones and Mitchell. We suggest that this stress field is relaxed by the 
mechanism which they have described, namely, the formation of positive 
prismatic dislocations which glide away over the surfaces of twelve 
slip cylinders with (110) directions radiating from the central particle 
as axes. The decoration of these prismatic dislocations with gold during 
the subsequent exposure of the crystals at room temperature and their 
heating to 140 to 160°C gives rise to the systems of regularly spaced fine 
precipitates which are shown in fig. 2, Pl. 9. The structure of these 
precipitates cannot be resolved with the optical microscope. 

Harvey and Mitchell (1958) have recently concluded that the systems of 
prismatic dislocations, formed during the cooling of crystals of sodium 
chloride containing aurous chloride, which have been described by Barber 
et al. (1957) and are illustrated in fig. 3 of their paper have the same origin. 
They have demonstrated that the dark central particles such as that shown 
in fig. 3(b) of the paper are specks of gold. The positive character of the 
prismatic dislocations has been established by observing their inclinations. 
in known stress fields (see fig. 7 in the paper by Jones and Mitchell 1958). 

Parasnis and Mitchell (1958) have also observed systems of prismatic 
dislocations extending along all twelve ¢110) directions radiating from 
particles of photolytic silver formed during the exposure of crystals of silver 
chloride containing 0-1 mol, of cuprous chloride. From these observa- 
tions they have concluded that space for the separation of silver is made 
available at the interface between a silver particle and the silver halide by 
the formation of positive prismatic dislocations which glide away from the 
interface along slip cylinders with the twelve (110) directions as axes. 

There can be little doubt that these groups of independent observations 
provide evidence for a previously unrecognized phenomenon in precipita- 
tion reactions in the solid state. When the rate of growth of a precipitate 
particle exceeds the rate at which the material of the matrix can be removed 
from the interface by the condensation of vacant lattice sites, a compressive 
stress field is established which is relaxed by the formation of positive 
prismatic dislocations, which glide away from the interface. We believe 
that this may prove to be of importance for the understanding of precipita - 
tion hardening because we have observed that the systems of closely 
spaced prismatic dislocations which result from mutual interactions when 
there is a high density of precipitate particles present an effective barrier 
to the glide motion of dislocations. The observations show that the 
strain associated with a particle of a precipitate can be transmitted to 
distances many time the radius of the particle by the mutual repulsion of 
the prismatic dislocations. It seems that the possibility that some 
features of the observations of Preston and Guinier on the diffraction of” 
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x-rays during the early stages of precipitation hardening in certain alloys 
could be interpreted in terms of diffraction by arrays of positive prismatic 
dislocations or by precipitates nucleated upon them should be considered, 


§ 4. SUMMARY AND CONCLUSION 


Dislocations may be decorated by the precipitation of gold throughout 
the volumes of sheet crystals of silver chloride with thicknesses between 
200 and 400. The essential features of the new experimental method are 
as follows. After separation from the glass plates, the crystals are cut 
into sections and annealed for 8 to 12 hours at a temperature between 
350 and 400°c in chlorine at a pressure between 200 and 300mm Hg 
During this anneal the areas around the cut edges of the crystals recrystal - 
lize and the remaining central area polygonizes. The crystals are then 
heated in a sealed tube in chlorine at a pressure of 15 to 30mm Hg and at a 
temperature of 350°c together with a piece of gold wire and 5x 10->gm 
of cupric chloride for 4 hours. They may then be strained at room 
temperature. They are exposed to light from a mercury vapour lamp 
to produce nuclei along the dislocation lines and heated at 140 to 160°c for 
1 to 15 minutes to increase the size of these nuclei until they are visible 
in the microscope. 

A new class of precipitation phenomena is described and discussed. 
Four independent groups of observations show that the compressive 
stress field around an inclusion or a growing precipitate particle may be 
relaxed by the generation of positive prismatic dislocations which glide 
away from the interface between the particle and the matrix along slip 
cylinders with the twelve (110) directions passing through the centre of 
the particle as axes. 
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ABSTRACT 


Measurements are reported of the electrical resistivity of pure poly- 
crystalline specimens of bismuth and antimony over the range 2 to 300°K 
and of the thermal conductivity up to 150°x. Some measurements in a 
magnetic field have also been made which allow separation of the thermal 
conductivity into electronic and lattice components. Below 50°xK, it 
appears that the heat conductivity of bismuth is nearly all due to lattice 
waves; free electrons seem to contribute little to the conduction or to the 
scattering. In antimony, however, the lattice conductivity is smaller than 
might be expected at low temperatures suggesting significant scattering by 
free electrons. The electronic term, A,, in antimony is remarkably similar 
in its temperature dependence to that in a monovalent metal. 

The electrical resistivity of both elements varies as about 7” for 7 < 6/10. 
For bismuth the approximate proportionality of resistivity with temperature 
observed at normal temperatures extends to rather lower temperatures than is 
usual in a metal. 


$1. InTRODUCTION 


BisMuTH and antimony are both elements with five valence electrons: 
which are capable of filling a so-called Jones zone exactly, but due to the 
smallness of the energy discontinuity at the zone boundary, a small 
number of electrons overlaps into the next band, leaving holes in the 
almost filled band (for example as discussed by Wilson 1953, pp. 92, 169 
or Mott and Jones 1936, p. 302). Measurements of such properties as 
magneto-resistance, Hall effect, diamagnetism, anomalous skin effect 
and de Haas—van Alphen effect enable estimates of the number of free 
electrons per atom to be made. Unfortunately these measurements 
(see Wilson 1953, p. 230 and Shoenberg 1952) do not all indicate similar 
values for this number, estimates for bismuth varying from as low as 
about 10°° to perhaps 10-? electron per atom. It appears that much 
experimental evidence suggests about 10~° electron per atom in bismuth 
and perhaps 10°? in antimony. At ordinary temperatures the electrical 
resistivities of these elements are only an order of magnitude higher than 
for many good metallic conductors; it thus seems that the electron 


+ Communicated by the Authors. 

{ The use of the name Jones zone for certain energy zones which are not 
strictly zones in the sense defined by Brillouin, has been recently discussed in 
detail in a review by Reitz (1955). 
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mean free paths in these semi-metals must be very much longer than in 
normal metals. Sondheimer (1952) has pointed out that a long mean 
free path should be a consequence of the small number of free electrons 
in the almost empty conduction band ; these electrons may only interact 
with lattice vibrations of long wavelength and such interactions will be 
relatively unimportant except at low temperatures. This might also be 
expected to lead to a linear variation of electrical resistance with temper- 
ature down to rather low temperatures, e.g. to about 10°K for bismuth. 

The small number of available free electrons may be affected appreciably 
by the presence of certain impurities, but fortunately antimony and 
bismuth are now readily available in a state of high purity and may be 
purified in the laboratory by zone melting techniques. Their electrical 
and thermal conductivities at low temperatures have been investigated 
experimentally for many years but most determinations have been 
concerned with the effect of magnetic fields and crystal anisotropy on 
these properties in rather restricted temperature intervals. 

We have measured the electrical resistance of our samples over the 
relatively wide temperature range extending from 2°K to 300°K and the 
thermal resistance from 2°K to about 150°k. Polycrystalline specimens 
of high purity and fairly large grain size were used. As well as estab- 
lishing the temperature dependence of the electrical resistance and of the 
thermal conductivity, we have tried to separate the latter into electronic 
and lattice components. 

References to previous measurements of thermal conductivity by 
Gehlhoff, Rausch, Rosenberg, Eucken, Reddemann, Rodine, de Haas, 
Gruneisen, Capel and others may be found in the compilation of Powell 
and Blanpied (1954) and the review by Klemens (1956). Data on electrical 
resistance is discussed by Wilson (1953), and in the recent review by 
Gerritsen (1956). 

A preliminary report of part of this work was presented (White and 
Woods 1957 a) at the Fifth International Conference on Low Temperature 
Physics and Chemistry (Madison, Wisconsin, U.S.A.) and thermal 
conductivity results for two of the less pure bismuth samples (Bi | and 
Bi 2) were reported in an earlier paper (White and Woods 1955). 


§ 2. EXPERIMENTAL METHOD 


The cryostat and method of measurement have already been described. 
(White and Woods 1955), helium-filled gas thermometers being used 
for temperature measurement and a galvanometer-amplifier for observing 
the electrical resistance. 

Table 1 lists the following characteristic parameters of the specimens 
used in this work: approximate diameter, D; residual electrical 
resistivity, p,; electrical resistivity at room temperature, py95 ; and the 
value assumed for the Debye characteristic temperature, 4p (see Keesom 
and Pearlman 1956). 
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The specimens were rods about 6 cm long and connections were made 
by soldering to them—soft solder for antimony and Woods’ metal for 
the bismuth specimens—except Bi 2 which was cast with copper potential 
leads in it. Their purity and general history are as follows 


Table 1. Data for Specimens 


D Po P295 Y 

(mm) (u Qem) | (wu Qem) | (°K) 
Bi 1 2 104 136 120 
Bi 2 3 5-9 120 120 
Bi 3 ~3-5 2:0, 118 + 120 
Bid |5x5x2:5| 1-70 118+ | 120 
Bi 5 6x6 2-4 118 + 120 
Sb 1 4:3 2-5 0-057 47-7 200 
Big Neale 0-054 41-3 | 200 
Sb 2aJ ss 
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+ Because of the irregular dimensions of these specimens, this room tem- 
perature resistivity value was assumed (cf. Gerritsen 1956); the geometrical 
form factor of the specimen(s) was then calculated using this value together 
with the measured resistance of the specimen : ie. //d =R/p. 


Bi 1 was cast and cooled quickly in a glass tube ; it contained columnar 
crystals penetrating to the centre of the rod, 16 to 18 crystals being 
exposed on the circular section. Its original purity of 99-97°, when 
supplied by Mining and Chemical Products (London) was no doubt 
affected by having been previously handled and cast before this final 
casting. 

Bi 2 was cast in a brass former from granular bismuth of about 99-99%, 
purity supplied by the General Chemical Division of Allied Chemical 
and Dye Corporation. Due to rather slower cooling, the columnar 
crystals appeared larger, there being about six crystals per circular 
section. 

Bi 3 was zone-refined by Mr. Gaston Fischer of this laboratory starting 
with 99-999°, bismuth from Varlacoid Chemical Company of New York. 
It was subsequently annealed in a larger tube for several days at a 
temperature just below the melting point. The crystals were about 
1 cm long and had the lateral dimensions of the rod; these dimensions 
were somewhat irregular as it melted temporarily during annealing. 

Bi 4 and Bi5 were both cut from a zone-refined bar prepared by Mr. 
Fischer from Varlacoid bismuth. Bi4 had a triangular cross section of 
sides 5, 5 and 2-5 mm and contained about three large crystals. Bi 5, of 
square section, consisted of crystals 1 to 2 cm long and from 2 to 4mm 
wide, slightly inclined to the axis of the rod. 

Sb 1 was sawn from a lump of ‘extra high purity’ grade antimony 
of the Bradley Mining Company (San Francisco); it was not annealed. 
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After measurements were made it was annealed but unfortunately broke 
before being remounted in the cryostat. However a second antimony 
specimen, Sb 2, was made by Mr. Fischer in this laboratory. It was 
prepared by zone-refining high-purity grade antimony from Bradley 
Mining Company and later annealing the specimen at about 600°c for a 
week. The crystals in Sb2 varied from about 2 to 5mm in width. 
The low electrical resistance ratio, po/py9,=1:3 x 10-3, is a confirmation 
of the high chemical purity of these specimens. 

Sb 2a is identical with Sb 2, having been merely remounted in the 
cryostat for the magneto-resistance studies. 


§ 3. RESULTS 
3.1. Thermal Conductivity 


Our results for bismuth, shown in fig. 1 together with values from 
Shalytt (1944), indicate that at temperatures below 3°K the heat 
conductivity, A, increases with temperature as 7’ where n=2-7--0-2; 
at higher temperatures, from about 10 to 50°x, \ falls approximately 
inversely as the temperature. As we shall discuss below, this behaviour 
seems quite typical of a pure crystalline solid in which lattice waves 
conduct practically all the heat unimpeded by the presence of any free 
electrons. Only at the highest temperatures is there a departure 
suggesting that electrons may be making a contribution to the con- 
ductivity. The change in thermal conductivity with magnetic field 
tends to confirm this. At liquid oxygen temperatures a field of about 
104 oe reduces the heat conductivity by 15 to 20% but has no noticeable 
effect at liquid helium temperatures, a behaviour observed also by 
Shalyt. The increase in electrical resistance due to a field of 104 0e is 
about 30-fold at 90°K, so that a field of this magnitude should be quite 
sufficient to reduce the electronic thermal conductivity, A,, to insigni- 
ficance compared with the lattice conductivity, 4,, provided that the 
Lorenz ratio , L,=pA,/T7’, is not seriously affected by the field. 

The results of our measurements on antimony are shown in fig. 2 
together with those of Rosenberg (1955) and of Gehlhoff and Neumeier 
(1913). These latter results from 80 to 300°K seem representative of 
much of the earlier work in this temperature range (see for example 
curves in the compilation of Powell and Blanpied 1954). In contrast 
to the behaviour in bismuth, our antimony results show that at very low 
temperatures }~aT' ; at least for Sb 2, the constant, a, has a value close to 
L|po= 2°45 X 10-8/p, W cm-1deg~?, which indicates that the electronic 
component, A, is dominant. For the unannealed rod Sb 1, the measured 
Lorenz ratio or the proportionality constant, a, is about 10% lower than 
the theoretical value. It seems most unlikely that this measurement is 
correct because it is well established for pure metals at very low temper- 
atures that pA,/7’—=L)o2 2°45 x 10-8W Q deg? and the presence of a lattice 


+ Obtained on a Hilger single crystal of high purity, which had a cross section 
of about 0-1 em2; the residual resistance ratio py/po; was about 20 x 10~, not 
very different from that for our samples 3, 4 and 5, 
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component only serves to increase the apparent value of Lo. As specimen 
Sb 1 was subsequently annealed and broken, another determination 
cannot be made but it is probable that the measurement of p, was in 
error. The specimen was short and of fairly large cross section so that 
an electrical resistance of about one micro-ohm had to be measured ; for 
this it was necessary to use the galvanometer amplifier at its maximum 
sensitivity and a small calibration error may have occurred. 


Fig. 1 


THERMAL CONDUCTIVITY (W cm-' deg -') 


| 2 = 10 20 50 100 
TEMPERATURE (°K) 


Thermal conductivity of bismuth ; experimental points are omitted for clarity 
but curves have been drawn through the points (cf. fig. 2). X represents 
values obtained by Shalyt (1944). 


As in the case of bismuth, the small effective number of free electrons 
in antimony is associated with a large magnetic field dependence of the 
electron transport properties. Since the electron component of thermal 
conductivity will be affected by magnetic fields while the lattice conduct- 
ivity should be field insensitive, this provides a method of separating 
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the ternis, A, and , The increase in the electrical resistance of Sb 2a at 
temperatures of 4:2°K, 11-4°x, 17-2°K, 24°K and 34:2°x when a field of 
5 koe was applied was 200-fold, 300-fold, 60-fold, 11-fold and 3-fold 
respectively. At temperatures below 20°xk, a field of 2 or 3 koe 
decreased the thermal conductivity considerably but larger fields, up to 
10* oe, had no measurably greater effect so we have assumed that this 


‘saturation ’ value of A, called Aj, is the lattice conductivity, A, (in fig. 7). 


Fig. 2 


THERMAL CONDUCTIVITY (W cm-! deg-') 


| ce) 10 50 100 500 
TEMPERATURE (°K) 


Thermal conductivity of antimony. 

1 ——— Rosenberg (1955) 

2 — -— Gehlhoff and Neumeier (1913) 
5) 


At temperatures between 20 and 40°K a field of about 10* oe was used to 
obtain a value for A,=A,. At 58 and 91°K the method used by Rausch 
(1947, see also review by Olsen and Rosenberg 1953) has been followed 
to obtain the lattice conductivity. As shown in fig. 3, A and p at a 
given temperature are measured in a series of magnetic fields ; then it is 
assumed that the Lorenz ratio, L,=pA,/7’, is field-independent at that 
temperature and if 
A=),+A, 
==1,7') [Mai r ‘ 
a graph of A vs 7/p should yield a straight line with intercept A, at Tio=0, 
Our values for \, as well as for the total thermal conductivity, Ne Nee x 
and for \,=A—A, of specimen Sb 2a are shown in fig. 7. 
2B2 
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3.2. Hlectrical Resistivity 

‘The rather impure specimen, Bi 1, had a very high residual resistance 
(see table | and White and Woods 1955), showed strong departure from 
Matthiessen’s additivity rule and a marked minimum in the electrical 
resistance at about 18°k. As discussed, for example, by Wilson (1953, 
p. 228), the band structure of bismuth is very sensitive to the presence of 
any impurities which increase or decrease the number of effective electrons 
or holes and these effects may be temperature sensitive giving a form of 
semiconductor behaviour in some cases. 


Fig. 3 
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A graph of thermal conductivity, A, as a function of 7'/p in magnetic 
fields of 0 to 14 koe. 


For the other bismuth and antimony specimens the electrical resistivity 
p, appears sensibly constant below about 3°K, and values of the residual 
resistivity, py, assumed to be due to scattering by static imperfections, 
are given in table 1. Because of an uncertainty of a few per cent in their 
geometry the resistivities of specimens Bi 3, 4, and 5 have been normalized 
to give a room temperature value, py,,=118 x 10° Q cm (see review by 
Gerritsen 1956). The values for the ideal resistivity, p,=p—p,, due to 
scattering by thermal vibrations, agree fairly well for different bismuth 
specimens, the maximum deviation being a few per cent between Bi 2 
and Bi 3 in the region from 30 to about 150°K. This presumably reflects 
the different average crystal orientations in these two specimens. Mean 
values of p, are given in table 2 and these accord quite well with mean 
values of the reduced electrical resistance, R/R,,,, for different single 
crystals of bismuth investigated by Schubnikow and de Haas (1930). 
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Table 2. Smoothed Values of ‘Ideal’ Electrical Resistivity in 
micro-ohm cm 


T°K Bi Sb Pein ©) 3h Sb 

6  44(120°K) 25-9 (200°x) 50. «19 3-2 
295 116 41-3 40) 15 2:1, 
273 105 37-6 30 10 1-2 
250 96 34-0 20. 5:8 0-4, 
200 74 25:9 15 3-4 Oars 
15055 17-9 10 1:3, 0-06, 
100.37 10-0 5 0-2, ~0-009 

75 29 6-5 

Fig. 4 
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Tdeal electrical resistivity, ppp —pp, of bismuth and antimony. Experimental 
points are omitted for clarity. 
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The form of p, as a function of temperature is illustrated in fig. 4. 
The rather high value shown in table | for the room temperature resistivity 
of Sb 1 (cf. review by Gerritsen 1956) may be due to uncertainty in the 
shape of the specimen, so the values quoted for antimony in table 2 were — 
taken from the results for both Sb 1 and Sb 2, the Sb 1 values having — 
first been normalized to give pog,=41-3 x 10~° $2 cm. 


§ 4. Discussion 
4.1. Thermal Resistivity 
General. “xperiment and theory on the thermal conductivity of — 
solids have not always been in good agreement during the past few years, 
but patterns of behaviour have emerged for crystalline insulators and 
for pure metallic elements which are fairly consistent. Reviews by 
Olsen and Rosenberg (1953), Berman (1953) and Klemens (1956) discuss 


the situation in detail. Briefly, we may expect the following behaviour : 
(a) For pure metallic elements 
A=A,+A, where A, >A, 
and A= W=W,=W t+ VW; : 
where the impurity resistivity, Wy=py/ll’. Also for 7 >6), the ideal 
resistivity, W,, becomes sensibly constant and equal to W,,, the total 
thermal resistivity at high temperatures. It has been shown (MacDonald 
et al. 1956, White and Woods 1957 b, White 1956) that to a good approx- 
imation the ratio W,/W ,, is well represented by the function : 
rOfD ade 
fot )(l ee) 
for most metallic elements. For 7'’<@, this expression reduces to 
W |W 214-4 (1/6p)?. 
(0) For crystalline insulators A=A, and 
1/A,= W=W,tWpt+WptWz, 
where the resistivity components arise from scattering of the lattice 
waves by the crystal boundaries (W,), by dislocations (Wp), by point 
imperfections (Wp), and from an anharmonic interaction between the 
lattice waves (W,). For a crystal with no imperfections the thermal 
conductivity is limited at low temperatures by boundary scattering and 
Aa T* while at high temperatures, where the lattice wave interactions 
become important, A=1/W,,«1/T. 

In metals a term, W,, must be added to W, and W, to account 
for scattering of the lattice waves by free electrons. This term, which 
varies as 7’ ” at low temperatures, is large enough to reduce A, to 
insignificance in comparison with A, if the metal is a pure element. In 
semiconductors and perhaps in semimetals (e.g. Bi, Sb), where the number 
of free electrons is small, this term is much less important. 

Bismuth.  Kesults on bismuth (fig. 1) suggest that boundary scattering 
of lattice waves is the main resistive mechanism at liquid. helium 


W,|W o=2(T/0p)? | 
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temperatures. The formula obtained by Casimir (1938) may be reduced 
to Ap=2-3DT+ for the case of bismuth where D is the effective diameter 
of the crystals. A comparison of effective diameters calculated from 


Fig. 5 
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Thermal conductivity of bismuth 
——— (Total) Bid 
— -— A, eale. from A,=Ly!'/p 
A, interpolated assuming W;=2(7'/6)*J 3(6/7’) 


x A, deduced from Shalyt (1944) 
o A, from Bi 5 
CJ X, from Griineisen et al. (1950) 


© A, from Reddemann (1934) 

+ The numerical factor in this equation is slightly different from that quoted 

previously (White and Woods 1955) as we are now using the value 6)—120°K 
suggested by recent specific heat data. 
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A (at 2°K) using this equation and the actual grain dimensions is given 
below. Agreement is far from perfect but in view of the uncertainty of 


the average grain dimension, it seems to confirm that grain boundaries — 
are the main cause of scattering of lattice waves at low temperatures. . 
a 


Table 3. Crystal Dimensions for Bismuth Samples. 


Sample 


D (grain size) cm | ~0-02| ~0-1 | 0-1-0:3 | 0-2 | 0-2-0-4 
D (from Avex) em | 0-02] 0-06 | 0-05 0-07 | 0-1, 


In fig. 5, the total observed conductivity of Bi5 is plotted again, 
together with the probable form of the lattice conductivity, A,, and 
possible form of the electron contribution, \,. The extrapolation of A, 
has been made with the aid of our measurements on Bi 5 in a magnetic 
field which gave \,~0-16 (91°K), 0-18 (78°K), 0-23 (60°K) and also those 
given by Reddemann (1934). 

Thus we find that A,o 7-11 which is not very different from the 
expected 1/7’ variation due to anharmonic coupling. If, in fact we 
approximate the data in fig. 5 by 


A,~14/T W cm deg, 


we may then compare this expression with theoretical ones deduced by 
Leibfried and Schlomann (1954), Dugdale and MacDonald (1955) and 
Kontorova (1957). These expressions, strictly speaking only valid for 
simple isotropic solids, are based on extensions of a linear chain model 
and are essentially equivalent to one another. If we assume the average 
interatomic distance in bismuth to be about 3-3 A and the Debye temper- 
ature to be 120°k, then the Leibfried—Schl6mann formula reduces to 
Ap=55/y*7'. It seems misleading to give even an approximate value for 
the Griineisen parameter, y, but we may note that if y=2 (as it does 
approximately for many cubic solids) then remarkably good agreement 
with the experimental expression is obtained (cf. our remarks on Ge and 
Si (White and Woods 1956) and Te (Fischer et al. 1957). 

Between 10 and 5°K the conductivity rises much more rapidly with 
decreasing temperature than a 1/7’ law would predict. Peierls (see for 
example reviews by Berman 1953 and Klemens 1956) has ascribed such 
a rise in the lattice conductivity to a rapidly decreasing probability of the 
Umklapp processes which gives rise to the thermal resistivity, W,. If 
we assume that Wp, Wp and W, are all zero so that 1/A=W,+Wy, we 
may extrapolate Ay = 1/Wy, as shown in fig. 5, down to about 4°x 
with some confidence. If we now write A~Cvl/3 where C is the specific 
heat per unit volume and v is the velocity of the phonons, (say 
1-5 x 10° cm/sec in Bi), then the mean free paths, / (shown in fig. .6) 
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can be calculated from A or from A, The apparent exponential depen- 
dence of [(Ay) upon 6/7 (fig. 6) which gives some confirmation of Peierls’ 
theory of Umklapp processes, has also been observed for solid helium, 
diamond and sapphire (see review by Berman 1953). 


L(Xy) 212x107 exp (8/2.8T) 
lo-! 
E Wes 
Tee 
oa 


@/T 


Mean free paths, J, of phonons calculated for Bi 5. 
Points © are calculated from Ay and points ® from A. 


Turning to the electronic component of heat conductivity, A,, in 
bismuth, we find less precise information. The low and high temper- 
ature extremities of the A, curve drawn in fig. 5, have been calculated 
from values of the residual electrical resistivity and the high temperature 
values of the ideal electrical resistivity (where p,oc 7’) assuming 


pA, |T'= 2-45 x 10-8 WQ deg”? 
It does seem that this Wiedemann—Franz—Lorenz law should be valid 
in regions where scattering is elastic, that is where the effective mean free 
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paths of the electrons are the same for both electrical and thermal trans- 
port processes. However at intermediate temperatures the behaviour 
of the ideal thermal resistivity, W, = 1/A,, for a semimetal like bismuth, 
seems unpredictable at present. In bismuth below 20°K the heat 
conducted by the electrons is so small compared to that conducted by 


the lattice that direct experimental measurement using a magnetic — 


field only gives a crude value. Above 50°K we might expect values of 
\,=A—A, to be more reliable and they fall rather below the values 
calculated from the W-F-L law except at room temperature where 
Reddemann’s measurements give values of about 0-06 W cm~! deg™?. 


Fig. 7 
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Thermal conductivity of antimony 
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——-——- A,=A—A, (Sb 2a) 


Antimony. ‘The measured total heat conductivity, A, for Sb 2 (and 
Sb 2a) is shown in fig. 7 together with values of the lattice conductivity 
measured with the aid of a magnetic field, and a dashed curve representing 
\—A,=A,. The isolated experimental values for \, and A, at about 
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80° and 90°K are those obtained by Rausch (1947) for a single crystal 
measured in two directions in a plane perpendicular to the axis of trigonal 
symmetry. ‘These two directions, which Rausch denoted by §, and §, 
and which gave similar results, are the [110] and [112] directions in the 
smallest rhombohedral unit cell. 

Consider first the lattice conductivity, A,, for antimony: At high 
temperatures A, = 34/7’ W cm~! deg~! which we may compare with a 
value of Ay = 110/y?7' deduced from the Leibfried-Schlbmann formula 
using 3-1 A for the average interatomic distance. We may again remark 
that a value of y in the vicinity of 2 gives close agreement but this should 
not be taken too seriously. The accepted values of cubic expansion 
coeflicient, «, and compressibility, y, for polycrystalline material lead to 
y=u/xU = | (ef. Mott and Jones 1936, p. 318) for both antimony and 
bismuth. However Barron (1956) has pointed out that for anisotropic 
materials in which 7 is small the variation in y among individual normal 
modes can not be neglected; even in the limit where y=0 there will, in 
general, still be cubic terms in the crystal potential causing coupling of the 
normal modes of vibration. 

At lower temperatures the behaviour of the lattice conductivity is 
rather more puzzling. In this very pure, well annealed sample, Sb 2, 
A, varies as about 7° below 4°k and the maximum mean free path of 
the lattice waves appears to be about 0-1 mm which is very much less 
than the grain size. Although we have no measurements in a magnetic 
field for the specimen Sb 1, comparison of Sb | and Sb 2 indicates that the 
lattice conductivity is much smaller still for Sb 1 at liquid helium temper- 
atures but a high density of dislocations (e.g. N ~101' dislocations per cm? 
according to the formula of Klemens 1956) would explain the low lattice 
conductivity in this case. In Sb 2 not only are we faced with a small 
lattice conductivity below 4 or 5°K but also with a rather flat maximum 
at ~10°K and no very rapid variation of A, with temperature just above 
10°K, i.e. between 6/20 and 6/10. This latter behaviour is similar to that 
in Ge, Si and Te but contrasts with that in bismuth or in solid helium, 
sapphire and diamond. It is tempting to attribute this flattening to the 
presence of other isotopes which contribute a point defect scattering 
term, W,, to the thermal resistivity, a term sufficient to mask the sharp 
maximum such as is observed in bismuth. Bismuth consists of one 
stable isotope whereas antimony consists of two: Mass 121 (56%) and 
mass 123 (44%). However, a calculation of the isotope scattering term, 
W >, using an equation due to Klemens (1957) yields Wp >! = Ap=1000/7' 
watts per cm deg. It appears, therefore, that the presence of more than 
one isotopic species cannot be responsible for limiting the maximum. 
Perhaps there are sufficient free electrons present in antimony to modify 
the lattice conductivity at low temperatures although in bismuth there 
is no evidence of such a phonon-—electron scattering process. 

Consider now the electronic thermal conductivity of antimony as 
exemplified by specimen Sb 2 or Sb 2a. It appears similar to the thermal 
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conduetivity of a metallic element in which there is about one free electron 
per atom. Assuming the validity of Matthiessen’s rule for thermal 
resistivity we may write 


1/A, = WV = Wot W ,. 


Now at 2°x, W,< W, and from values of the electronic thermal con- 
ductivity and the electrical resistivity, py, in this temperature region we 
obtain po/W)T' 2-3, 10-8 W Qdeg for Sb2. Assuming that the 
factor, W,7', remains constant as the temperature is increased we have 
calculated values for the ideal thermal resistivity, W,;. As shown in fig. 8, 
W, appears to be approaching a constant value of W,,—7 cm deg W =! 
at high temperatures, which compares with a value of about 6 cm deg W 
calculated from the electrical resisitivity near room temperature, using 
the Wiedemann—Franz—Lorenz law. 


W, (cm deg w-!) 


5 10 50 100 500 
TEMPERATURE (°K) 


Values of ‘ ideal ’ electronic thermal resistivity, 


W;, deduced for antimony. The broken curve represents 
Wi =2W (7/6) 5(0/7’) assuming 6=200°K, and W..~7-5 cm deg W-1, 


At low temperatures, W,(7') is not very different from the function 
mentioned earlier :- 


oz a? da 


Mesa | ee 


=2W (7/0)? J5(6/7'). 
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This semi-empirical function gives a useful representation of data for the 
alkali metals and many divalent and transition elements but it seems a 
little surprising that it may be extended to cover, even in an approximate 
fashion, data for antimony. 


4.2. Hlectrical Resistivity 
General. For a simple metal in which the charge carriers are free 
electrons in a single band the Gritneisen—Bloch relation predicts 


Ts pe xv dx 
o | 


oe e°—1)(l1—e-+) 


pix = (T°/0°) J5(0/T) 
which for 7’< 6 reduces to 
Asoc). 


If the characteristic temperature, @, is identified with the Debye tem- 
perature, Op+, good agreement with experiment is obtained in the case 
of sodium. For other monovalent metals the 7° dependence is usually 
observed at very low temperatures, but a value of 6 somewhat less than 
8p) may be required to give good agreement with experiment for the ratio 
of low and high temperature values. 

For the transition elements in which the charge carriers in two bands 
(s and d bands) have very different effective masses, Wilson (1938) 
has shown that, provided interband transitions are not forbidden by 
momentum conservation requirements, there should be a large additional 
resistivity, 

oe 2 Os 


pea (POP | eet ew SUPP TOP) 


Kerei< 0) psqocl”. 


In elements like antimony and bismuth, requirements of momentum 
conservation have a further effect (Sondheimer 1952). Only those 
lattice waves with wave number, q, smaller in magnitude than 2|k| may 
interact with electrons having wave vector, k, and for the electrons at 
the surface of the Fermi distribution k will be small in bismuth and 
antimony. Hence at normal temperatures the small number of free 
electrons may only interact with a relatively small group of phonons ; 
thereby they have a large mean free path and also the range over which 
pxT is extended to quite low temperatures. 

Bismuth and Antimony. In both these elements, our data for the 
polycrystalline specimens (fig. 4) indicate that for 7’<6/10, paT®”. 
For bismuth, however, the transition from an approximately linear 
dependence at higher temperatures to this 7?’ dependence is much more 
sudden than for antimony. For neither element do the results appear 
to follow a Griineisen—Bloch function, but for antimony p,; may be fitted 


+6, as determined from specific heat data obtained in the range of tem- 
perature between 7'~ @p/2 and T'=@p. 
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tolerably well by the (7'/@)8 J,(6/7) function assuming 0=0)—200°K. 

Thus the ratio, p,(4/10)/p,(@), has values of 0-005 for the Grimeisen— 
* Bloch function, 0-015 for the Wilson function, while experimentally it 
has the values 0-015 in Sb, and 0-046 in Bi (assuming 6=6)—120°K for 
bismuth). This at least suggests that interband transitions are important 
in antimony and perhaps also in bismuth but are modified in bismuth by 
the smaller number of free electrons which cause an extension of the 
temperature region over which elastic scattering occurs, i.e. over which 
pa Ls 

§ 5. CONCLUSIONS 

(a) For pure bismuth, the heat transported by the lattice waves is very 
much greater than by the electrons (A, >A,) at temperatures below 50°K, 
and the lattice thermal conductivity is limited by the presence of grain 
boundaries and Umklapp processes in much the same way as it is limited 
in a pure dielectric crystal. No useful information about the minor 
component, ,, can be deduced as yet for the interesting temperature 
range below 50°K. 

(b) For pure antimony the two components of conductivity are com- 
parable except at very low or high temperatures. The electronic com- 
ponent, A,, behaves very like that in a simple monovalent metal. The 
lattice component, A,, at temperatures above 20°K, appears to be limited 
by anharmonic interaction between lattice waves but at low temperatures 
is much smaller than we should expect for a crystalline insulator (or than 
is observed for bismuth) suggesting that free electrons may seriously 
scatter the lattice waves. 

(c) The electrical resistivity of both elements varies as about 7275 
below 6/10 and depends linearly on 7’ at high temperatures. For bismuth 
this linear region seems to extend to considerably lower temperatures than 
for antimony. Interband scattering may be important in determining 
the temperature dependence of the electrical resistivity. 
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ABSTRACT 


It is shown that the polarization of luminescence from a cubic crystal, 
when excited by plane polarized light, may be used to determine the sym-: 
metry axes of the luminescent centres. This information is most readily 
obtained from a specimen which is bounded by (100) faces and by using light 
travelling normally to the faces. Experiments on a diamond cut in this way 
are described. The results are for the green luminescence from centres 
induced by heat treatment after irradiation. These centres appear to have 
symmetry axes along (110) directions, 


§ 1. INTRODUCTION 


Many of the optical properties of cubic crystals are isotropic, and it is 
not then possible to obtain any information about the symmetry of the 
centres which give rise to these properties. The polarization of light 
emitted by centres excited by plane polarized light is, however, anisotropic. 
This effect was first reported by Feofilov (1953, 1956) who investigated 
various centres in CaF,, LiF, and NaF. In the present note, measure- 
ments on diamond will be described, the method used being somewhat 
different. The optimum crystal directions were chosen on the basis 
of the theory of the effect, the relevant results of which are summarized, 
and from the results the axes of the centres are determined. 

There are two main luminescent centres in diamonds, the emission 
being in the blue and the green respectively. The latter centre can 
be induced in type I diamonds by heat treatment after electron, gamma 
or pile irradiation. Clark ef al. (1956) found that the absorption spectrum 
associated with this centre, when measured at liquid air temperature, 
consisted of a sharp line at 5035 A with a fine structure of other lines 
superimposed on a broad band on the high energy side of the line. The 
emission spectrum is similar (Dyer and Matthews 1958) except that the 
structure is on the low energy side of the main emission line at 5032 4. 

In this paper we report measurements on the green emitting centre 
when stimulated with 3650 4 and 4358 A mercury radiation. Experiments 
are continuing both on this and on the blue centre under varying conditions 
and these will be reported later with a fuller theoretical discussion. 


+ Communicated by the Authors. 
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§ 2. EXPERIMENTAL 

The apparatus is shown diagrammatically in fig. 1. When using 
3650 A excitation, light from a high pressure mercury lamp passed 
through a black glass filter (Chance—OX1) which cut out the visible 
light. In 4358 4 excitation the light from the lamp was dispersed by a 
f/8 glass prism Littrow spectrograph, the line selected by a slit and 
passed through three OB10 filters to remove a small amount of scattered 
green light. A Glan-Thomson prism P was used to polarize the exciting 
light which was then focused onto the diamond. The green emission 


Fig. | 
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Optical system. (S, source; P, Glan-Thomson polarizer; F, and F,, filters ; 
L, lens; D, diamond; R, rotating polaroid; and Ph.X, photo- 
multiplier.) 


Deflection (d) 


Time (angle of polaroid, /3’) 


Photomultiplier signal as a function of time. 


from the diamond was concentrated by the lens L and passed through 
two non-luminescent, green transmitting filters (OY4, OGrl). The 
emitted light was analysed by a continuously rotating polaroid and 
detected by a 14-stage photomultiplier. From the photomultiplier the 
output was fed either to a galvanometer or to a Brown recorder. . 
The light intensity recorded as a function of time (which is proportional 
to the angle of rotation of the analyser) was of the form shown in fig. 2. 


P.M, fi 
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In the experiments to be described the degree of polarization (p)— 
defined as the amplitude (a) of the oscillation divided by the mean value 
(/)—has been determined as a function of the angle of the electric” 
vector E of the exciting light relative to the crystal axes. 


§ 3. THEORY 


We begin by assuming the simplest forms of the symmetry of the 
centres in the crystal, namely that each centre has an axis of symmetry. 
If a typical direction of this axis is, say <A,, 4;, v;), there will be 23 other 
directions which are equivalent in a cubic crystal and there will be equal — 
numbers of centres having these 24 axes. In any particular centre with 
axis z, the energy levels may be singlets or doublets. For example, 
s and p, states would be singlets whereas the p,p, state would be a — 
doublet. Ifthe transition involved is an allowed electric dipole transition 
it may be of two kinds, one corresponding to what are usually called 7 — 
oscillators, the other to o oscillators. The first, of which the s—p, transi- 
tion is a typical example, has an intensity of absorption proportional — 
to cos? 6, where 6; is the angle between the axis of symmetry Z; and the 
E vector of the exciting light. The second has intensity proportional to 
sin? §, and the s—p,, p, transition is a typical case. The intensity of — 
absorption in the crystal is obtained by summing cos? 0,, or sin? 6; over 
all the equivalent directions. This sum is independent of the direction 
of E and the absorption is isotropic. 

The absorption of the centre in diamond was observed to be isotropic 
and we therefore conclude that the absorption is dipolar in character. 
Also, although it is not possible to find the number of centres accurately, 
rough estimates of the number of vacancies and interstitials produced 
in the irradiation indicate that the f-value of the transition is of the order 
of 10-1. We therefore rule out magnetic dipole and conclude that the 
transition is allowed electric dipole. 

If the light is emitted by the same 7 oscillator in which absorption 
took place the intensity of emitted light with E vector making an angle 
0,’ with axis Z,; is proportional to cos? 0; cos? 6,’._ If the direction cosines 
of the E vector of the exciting light are (/, m,n) the intensity emitted 
with E vector (l’, m’, n’) by the crystal is 


T=) (A, + mp, tr, PUA, mp, +n'y,)?, 6 6. (A) 


where J) depends on the experimental conditions and the properties of 
the centres (f-value and luminescent efficiency). This may be written 
I= fl y{ > PIP+ [YAP F I+ 2 (> WP—5 > PL ]}, - (2) 
lmn a linn Imn 
where the sums are over the three principal directions. For o oscillators 
the expression is very similar, 


T=Hy{1+ SPUP+( SAM PI +2 (SU P—5 YP}. (3) 
l 


lmn 4 mn linn 
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Asn? +ber?2+v2A,*?) is constant for the equivalent directions in a 
cubic crystal and its value characterizes the polarization of the lumin- 
escence. It varies from 0 for (100) centres to } for (111). For 
intermediate values there are an infinity of directions of the axes which 
give the same result. However, if the centre is of atomic size it is likely 
to have symmetry axes which are related to the atomic arrangement in 
the crystal, and are therefore in directions of high symmetry. For this 


reason centres with (110) axes are also studied for which YA?u?=H. 


({211) centres give the same value but are unlikely to occur in the 
diamond structure.) 

Many variations of the formulae are possible if more than two energy 
levels are involved in the process, and the system can make transitions 
between separate excited levels before fluorescing, or if the centres can 
change their symmetry axes when excited into high energy vibrational 
states. In this note, however, we shall simply test the applicability of 
the simple theory. 


§ 4. FoRMULAE FOR EXPERIMENTAL CONDITIONS 


In the experiment a rectangular block is used and the light enters 
and leaves the crystal perpendicular to two faces which may be parallel 
or at right angles. A direction normal to another face which is perpendi- 
cular to the light paths is chosen and the direction of the E vectors 
(l, m,n), (l’,m’,n’) are measured as angles f, 6’ from this direction. 
Equations (2) and (3) may then be written in the form 


I= M(B8)+2A(f) cos 26’4+2B(f8) sin 2p’. bie ew ver (4) 


In the experiment we cannot determine the intensity of light emitted 
having E at angle 8’—we can only use an analyser with its easy direction 
along f’ so that a fraction cos? 6, of emitted light having E at any angle 
6,, from the easy direction is also transmitted. Thus for any collection 
of dipole oscillators the deflection (d) cannot go to zero, the maximum 
value of p being 0-5 which would be obtained, for example, with a single 
dipole oscillator. 

It is convenient to express the deflection (d) corresponding to the 
photomultiplier current as 


d=M(f)-+a(f) cos 2 [B’—<(f)], - - - - + (©) 


where M(8), a(8) and <(8) are the mean, amplitude and phase of the 
oscillatory trace in fig. 2. The amplitude is related to the constants A 
and B by a(8)=—(A?+B?)!? and the phase by tan~? 2¢(8)=B/A. 

For rectangular crystals cut in a general direction the expressions for 
M, A and B are complicated and because the values of >A 7; are not 


very different for the (111) and (110) centres it is sometimes difficult 
to distinguish these two centres, 
202 
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The simplest expression for d and the maximum distinction between 
(100), {110) and (111) centres is obtained when the specimen has 
(100) faces. The angles 8 and f’ are then defined as shown in fig. 3 and 
from eqns. (2) and (3) we obtain the values shown in the table. Figures 
4, 5 and 6 show the variation of p with B for o and 7 oscillators under 
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Comparison of experimental results with the theoretical curves for the variation 
of p with B for ¢100), (110) and 111) centres, viewing perpendicular 
to the incident light. 


various conditions. It will be seen that the shapes of the curves are the 
same in the two cases, but that the values of p are less for the o oscillators 
owing to a constant being added to M. This depolarization arises 
because, when viewed along the z-axis the o oscillators emit circularly 
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polarized light, which leads to a constant contribution to the intensity 
transmitted by the rotating analyser. 


§ 5. Resutrs anp Discussion 


_ Dr. J. F. H. Custers of the Diamond Research Laboratory, Johannes- 
burg, kindly prepared for us a specimen in the form of a cube bounded, 
within 2°, by (100) faces. This crystal was given successive pile irradi- 
ation and heat treatment (800°c) until a sufficiently strong 5032 4 system 
was present. 

The crystal was birefringent along one direction but measurements were 
made for the other two straight through positions and the two right 
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Comparison of experimental results with the theoretical curves for the variation 
of p with B for (110) centres, viewing parallel to the incident light. 


angle positions not involving the birefringent direction. Because of the 
low conversion efficiency it was necessary to use a high pressure mercury 
lamp for the source. While the lamp was sufficiently stable to obtain 
reasonably accurate values of a/M for different values of 6, we were not 
able to obtain experimentally reliable values of a and M separately. 
Experiments were carried out using 4358, 4047 and 3650 A excitation. 
The results for 4358 and 3650 are shown in figs. 5and 6. Excitation with 
4047 A gave values of p equal, within the experimental error of about 
+0-005, to those obtained when using 4358 A—a result which was also 
found in another green emitting crystal whose faces were not cut in simple 


directions. 
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The experimental results in figs. 5 and 6 show a maximum value of p 
at P=45° in the straight through and at B=0 in the perpendicular con- 
ditions. In the former case this eliminates (100) centres for which the 
maximum would be expected at B=0. It is clear that for 4358 A excita- 
tion the results can best be described by using 7 oscillators in centres 
having (110) axes. 

We have discussed first the results for 4358 A excitation because they 
give the highest values of p. It is then necessary to consider how p 
could be reduced to account for the results using 3650 A excitation. The 
opposite approach is not possible because there do not seem to be any 
processes which would allow the degree of polarization from a given set 
of oscillators to be increased. 

Although the results for 3650 A excitation lie on the curves predicted 
for o oscillators about (110) axes there are two points which indicate 
that this is not the correct explanation. First, one of us (I. G. M.) has 
found that the fluorescence spectra obtained with the two excitation 
wavelengths are the same. Thus the emission must involve the same 
energy levels in the two cases, and we have concluded above that for 
4358 A excitation these energy levels are associated with « oscillators. 
Moreover, measurements on the other green emitting crystal which was 
referred to earlier, show that changing the excitation from 4358 A to 
3650 A can change the shape—as well as the scale—of the curves. We 
think, therefore, that some other process is present which gives a de- 
polarization when the excitation wavelength is 3650 A. 

The depolarization may arise because the centre is excited into a higher 
excited level than that from which emission takes place. Such a level 
would have to have a small transition probability for,excitation from the 
ground state since its absorption lines are not observed. In fact, the 
absorption at 3650 A, as measured by Clark et al. (1956), is small and it 
was only the high intensity of illumination available at 3650 A which 
made this a convenient source in preliminary experiments. The exci- 
tation might even be into the conduction band, as a result of which 
electrons could migrate from one centre to another and tend to destroy 
the polarization. 

If on the other hand the excitation occurs into the excited electronic 
state which has a minimum energy of 2-46 ev (5032 A), it does so to a 
state where the crystal has 0-93 ev of vibrational energy. Such excessive 
excitation might allow the centre to change its axes, thus giving less 
polarization. 

Concerning the atomic nature of the centre Clark et al. (1956) have 
suggested that it is formed when a product of damage—vacancy or inter- 
stitial carbon atom—migrates to some other imperfection which might 
be an impurity atom. The simplest double defect having (110) axes is 
an interstitial carbon atom anchored to an interstitial impurity. Com- 
binations such as vacancy /interstitial impurity, vacancy/substitutional 
impurity, or interstitial carbon/substitutional impurity would be expected 
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to exhibit strong (111) directional properties. The next simplest com- 
bination having (110) would be a vacancy anchored to a substitutional 
impurity in a next nearest neighbour position. 

Further experiments are in progress to enable measurements of the 
three parameters /, a and « to be made separately and to determine the 
cause of the depolarization with 3650 A excitation. The former seems 
to be entirely a problem of the stability of a source of the required 
intensity. Some gain should also be obtained by cooling the specimen to 
liquid nitrogen temperature, where a higher fluorescence efficiency is 
obtained (e.g. Dyer and Matthews 1957). 


§ 6. CONCLUSION 
The green emitting centre induced in a specially cut diamond by heat 


treatment after irradiation has been shown to have symmetry axes along 
«110). 
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ABSTRACT 


Observations of the intensity distribution of radio emission have shown 
a bright belt approximately normal to the galactic plane and passing 
through the anti-centre. Anexplanation is suggested in terms of the highly 
anisotropic radiation from relativistic electrons moving in the magnetic 
field of the local spiral arm. This suggestion avoids the supposition that the 
sun is in a special position in the galaxy. 


§ 1. [INTRODUCTION 


THE observed distribution of radio brightness over the sky has been 
interpreted in terms of symmetrical models of the galaxy by a number of 
authors including Hanbury Brown and Hazard (1953) and Baldwin 
(1955 a). The latter has proposed a model whose main components are 
(i) a spherical halo about the galaxy, (ii) a disc-like distribution of sources 
concentrated towards the galactic plane and its centre (Westerhout and 
Oort 1951), and (iii) an isotropic distribution of extra-galactic sources. 
Other less symmetrical features of the background radiation, which 
cannot be explained by the above model have been noted. The most 
important of these is a belt of emission from directions approximating to a 
great circle through the galactic poles and the plane at 1=0°, and 1=150°. 
It appears on a number of surveys, including among the most recent 
Baldwin (1955b), Piddington and Trent (1956), Blythe (1957). A 
diagram indicating the position of this feature is shown in fig. 1, and 
cross sections showing the angular variation of intensity at two points 
in the belt are given in fig. 2; the latter curves were obtained by Blythe 
(1957) on 38 Mc/s. The intensity of the belt exceeds the isotropic com- 
ponent due to extra-galactic sources and the feature is therefore likely 
to be of galactic origin ; it could be explained by an emitting disc, orientated 
normal to the galactic plane and containing the sun near its centre. If 
similar objects of the same brightness occurred elsewhere in the galaxy 
they would, however, be readily detectable individually, and would 
contribute the major part of the total galactic emission. No such objects 
have been observed, and it is the purpose of this paper to suggest an 
alternative explanation which avoids the assumption of any marked 
asymmetry of the galaxy. It is also an alternative to the proposal that 
there is a perturbing galaxy in collision with our own (Johnson 1957). 


y+ Communicated by M. Ryle, F.R.S. 
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Fig. 1 
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The broken line on the celestial sphere shows the belt’s approximate position 
relative to the sun S, the galactic centre C and the galactic plane. 


(i) ea (ii) mark the positions where cross sections have been plotted 
in fig. 2. 


Fig. 2 
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20 10 e) 
Declination 
Two cross sections of the belt at galactic latitudes (i) 28° and (ii) 70°. At this 
longitude a distribution in declination is effectively the distribution 
across the belt. 
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The plane of the apparent emitting disc is almost perpendicular to the 
mean direction of the local spiral arm in the neighbourhood of the sun. 
It is shown in § 2 that the intensity is not uniform around this disc; the 
variation shows a reasonable similarity to the cross sectional shape of the 
local arm as revealed by observations of the 21 cm line of neutral hydrogen. 
This similarity, therefore, suggests that the belt may be due to a source 
of emission situated within this local arm. 

The application of the synchrotron mechanism of radio emission 
(Shklovsky 1952) to relativistic electrons moving in a general magnetic 
field associated with the local spiral arm appears to present a possible 
explanation of the feature. It is shown in §3 that the belt might arise 
from the anisotropic radio emission which would be expected under such 
conditions. 


§ 2. THe Loca SprraAL ARM 


Observations of the 21cm emission at different latitudes in the plane 
passing through the centre and anti-centre, normal to the galactic plane 
(Heeschen and Lilley 1954) have shown that the cross section of the local 


Fig. 3 
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The distribution in galactic latitude of the maximum brightness temperature 
of the belt is compared to the cross section of the local arm as observed 
by Heeschen and Lilly (1954) at 21 cm. 


arm is not symmetrical; it shows a marked flattening which does not 
coincide with the galactic plane but is inclined at about 20° to it. These 
observations are plotted in polar coordinates in fig. 3, together with a 
similar plot of the variation of the maximum intensity of the belt as 
derived from the surveys at 38 Mc/s and 81-5Me/s. The latter sets of 


ss. 
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points do not correspond exactly to the plane of the 2lem observations, 
eo makes a small angle to the plane through the centre and 

The similarity between the curves suggests that the belt has its origin 
in the local arm and this possibility will now be examined. The existence 
within the arm of fast electrons having energies of at least 10®ev must 
be assumed}. The question of their origin is not yet answered with 
certainty. Following one suggestion (Fermi 1949 and 1954, Kaplan 1956), 
it will be assumed that they appear with the cosmic rays in regions where 
there are dense clouds in random turbulent motion. Near the sun, 
however, they are in the presence of an ordered magnetic field along the 
spiral arm for which there is quite independent evidence (Hiltner 1951, 
Davis and Greenstein 1951), the electrons would here spiral smoothly 
with a radius of curvature of the order of 2x 10-7pe. Their lifetime is 
limited by nuclear collisions to roughly 107 years but for this time they are 
effectively tied to the lines of force. This field, according to Fermi and 
Chandrasekhar (1953), is of the order of 6 x 10-* gauss and is regular in 
direction to within 5° over distances up to I1kpc. The regions where 
turbulent clouds exist and where the electrons are supposed to be accelerated 
are presumably further away than this. It will be supposed that such 
regions are confined to the inner parts of the galaxy entirely, in which 
case the electrons leave the inner turbulent region, follow the smooth 
field along the spiral arm and are eventually lost in the halo. 


§ 3. THE SuGGESTED MECHANISM 


The acceleration of fast electrons in regions of turbulent cloud motion 
will presumably give rise to an isotropic distribution of electron velocities. 
If such particles now diffuse outwards along the ordered magnetic field 
associated with the spiral arm, those having initial velocities nearly 
parallel to the field will move quickly into the outer parts of the galaxy, 
while those whose initial velocities are nearly perpendicular to the field 
will move comparatively slowly. At any point along the spiral arm 
there will, therefore, be many more electrons per unit volume having 
motions across the magnetic field than those with motions along it. 

Since the radio emission from a relativistic electron moving in a magnetic 
field is confined to directions very near the instantaneous velocity vector, 
the emission per unit volume will be stronger in directions perpendicular 
to the magnetic field than in directions along it. It can be shown then if L 
is the distance through which the electrons have drifted in the ordered 
field and 7 is their lifetime, then the intensity of the radio emission has a 
maximum at an angle cos-!(L/c7’) to the magnetic field, and that the 
radiation is restricted to an angular width of approximately L/c7' radians. 


esr pe eee eee 

+ This value is the average of the energies calculated as those chiefly respons- 
ible for the Cambridge Survey frequencies of 38 Mc/s and 81-5 Mc/s, assuming 
a magnetic field of 6 x 10-6 gauss, 
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If L=1kpe and 7'=107 years, the radiation is emitted at an angle of 
89° 59’ to the field with a half-width of 1’. A uniform field would therefore 
give rise to an extremely narrow disc of emission practically normal to the 
field direction. The actual configuration of the magnetic field of the 
spiral arm may be compared to a rough skein of wool, each strand represent - 
ing a line of force; the variation of emission from the whole cross section 
of the spiral arm will then correspond to the variation in the direction — 
of the magnetic field. Fermi and Chandrasekhar’s value for the latter 
is of the same order as the observed width of the belt. 

In order to account for the observed brightness of the belt, it must be 
supposed that the number density of fast electrons is about 20 times that 
in the halo (Baldwin 1955 a). 

If the halo is supplied by the drift of electrons from the inner part of the 
galaxy along the spiral arms, the ratio of the volume of the halo and the 
spiral arms would lead to a ratio of this order of magnitude. 

In addition to the radio emission from the local spiral arm, the mechanism 
suggested would also produce emission which would be detectable from 
regions in other spiral arms where the line of sight was perpendicular to 
the arm. The large increase in intensity observed at low latitudes within 
10° of the galactic centre and anti-centre may be explained in this way; 
the latter feature is particularly difficult to explain by any other mech- — 
anisms. 

It is also possible that some of the anomalies between the distribution 
of radio and optical brightness in certain external galaxies might be 
explained; a spiral galaxy would exhibit an additional component of 
radio emission in the form of a strip perpendicular to the major axis. 


. . 
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§ 4. Discussion 


It has been shown how this belt, an otherwise anomalous feature of 
galactic radio emission, may be due an isotropic distribution of the fast 
electron velocities near the sun, and how this distribution might arise. 
No special position need be assumed for the sun because a similar belt 
would be observed wherever the sun were placed within the spiral arm. 
The brightness and its variation with latitude are consistent with other 
data on the local arm. On the other hand, three major objections may be 
raised. (i) No similar anisotropy has been observed in the cosmic ray 
flux at the earth ; (ii) the direction of the local spiral arm does not coincide 
with the axis of the belt ; (iii) although synchrotron radiation is completely 
polarized, Thomson (1957) found less than 1° polarization in the belt. 

As for the first objection, one would expect cosmic rays observed at the 
earth to have been affected by the extremely local fields of the solar system. 
The emitting electrons are outside this region. These fields are sufficient 
to make the cosmic ray flux isotropic if they are of the order of 10-5 gauss 
(Richtmyer and Teller 1949). At energies of 10!” ev, however, it is inter- 
esting to note that a slight minimum has been observed as the earth 
faces inward along the spiral arm (Cranshaw and Galbraith 1954). 
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The second objection is that the plane of the belt is normal to longitude 
70°, whereas from 21em observations Westerhout (1957) has deduced 
that the Orion arm is directed towards longitude 65°. This discrepancy 
_and the finite width of the belt would readily be explained if the direction 
of the magnetic field were supposed to show variations of the order of 5° 
as shown in fig. 4. Since the radiation is emitted within a very small angle 
perpendicular to the direction of the local magnetic field, those regions 
which contribute to the received radiation lie on a curve such as C in 
fig. 4; the angular width of the belt would therefore be of the same order 
as the deviations in the direction of the field. 


Fig. 4 


Diagram to illustrate a possible type of local magnetic field. The source 
~ elements lie along the curve C, the locus of points at which the line of 
sight is normal to the lines of force. 


If this interpretation of the angular width of the belt is correct, the 
absence of appreciable polarization may be explained in terms of the 
Faraday rotation of the electric vector caused by the small longitudinal 
component of the field in the medium between the emitting element and 
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the earth. The emission from each element undergoes a rotation of its 
plane of polarization whose magnitude increases progressively with the 
distance of the element, since for simple fields of the type shown, the sign of 
the longitudinal component of magnetic field is constant. If the angle 
of rotation induced between the nearest and the most remote elements is 
appreciably greater than 7, the resulting polarization may be small. 
For an electron density of 0-1cm* within the spiral arm (Westerhout and 
Oort 1951) and the observed angular width of the belt of 5°, corresponding 
to an average longitudinal component of the field of 10~’ gauss, the angle 
of rotation for the frequency used by Thomson (160 Me/s) is of the order 
of 0-03 radians per parsec. 

It is therefore possible that for emitting elements extending over a 
depth of 100-150 parsec, the resultant polarization will be greatly reduced. 
If this interpretation is correct, observations at higher frequencies, or 
with sufficient angular resolving power to observe only part of the belt, 
might be successful in detecting an appreciable polarization. 
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ABSTRACT 


A lattice arrangement of G.M. counters has been used to measure the 
size spectrum of showers containing 3x 10% 3108 particles. In this 
region, & power law distribution is found, 

J(> N)=9-5 x 10-5(108/N)1:88+0-15 m- hr-! gter-}, 
An independent method of measuring the number of particles in the 
showers gave good agreement with these figures. 


§ 1. IyTROoDUCTION 


WE have used the apparatus described in I (Cranshaw and Galbraith 1957) 
to derive the size spectrum of air showers in the region 3 x 1083 x 108 
particles. The method consists in determining the size spectrum in 
terms of the number of units of the array which are triggered, and 
calculating with the aid of the density distribution measured in III 
(Cranshaw, Galbraith and Porter 1957) the corresponding size spectrum 
in numbers of particles. Finally an independent measurement of the 
mean number of particles in the shower triggering a given number of 
units was made, and this measurement is found to agree satisfactorily 
with the values calculated from the spectrum. This agreement gives 
confidence in the correctness of the assumptions in III. In particular, 
the number of particles in showers can be determined by density measure- 
ments at about 100-200m, and the use of the Nishimura and Kamata 
curve. 
§ 2. EXPERIMENTAL Data 

For the analysis, we wish to know the rate of showers which would 
discharge given numbers of units on an infinite lattice array with the 
spacing of our experimental array. For this purpose, we used the punched 
ecards showing the units struck on the array. In order to minimize edge 
effects, an area inside the array was marked out as shown in fig. 1. For 
>3 fold, > 5 fold or >8 fold showers, showers were accepted if the centroid 
of the units struck lay inside the smaller area. It is assumed that on the 
average the centroid of the units struck coincides with the core of the 


+ Communicated by the Authors. 
+ On leave from the University of Potchefstrom, South Africa. 


2 D 


378 T. E. Cranshaw ef al. on 


shower. For large showers, the whole of the array was accepted and 
visual estimates of the core of the shower were made. When a large 
part of the shower lies off the array, it is difficult to estimate how many 
units would be struck on an infinite array. Nevertheless, independent 
counts by different observers agreed well on the numbers of showers 
triggering > 20 and >40 units. Above 60 units, the uncertainty is large, 
but not larger than the statistical errors. The numbers of showers counted — 
range from 400 for >3 fold showers to 6 for > 60 fold showers. 


Fig. 1 


e GM units 


§ 3. DETERMINATION OF SIZE SPECTRUM 


To convert the size spectrum in numbers of units to the spectrum in 
number of particles, we first calculate the probabilities P(N, f) that a 
shower containing N particles will trigger f units. We calculate these 
quantities for the case of an infinite array with the spacing of our lattice, 
i.e. we ignore edge effects. 

For this calculation, we need to know the probabilities p, for each of 
the units to be triggered. If the position of the axis on the lattice is 
known, then the distances 7; to all the units are known, and we can use the 
density distribution function given in III to calculate the density at the | 
unit for a given N, and hence the probability of triggering. To obtain the 
value of N in III we have assumed the validity of the Nishimura and 
Kamata distribution for distances less than 50 metres, and therefore for 
the time being, we regard N as a nominal value, and not a measured value. 

We have used two methods to calculate P(V,f) from the p,. First, the 
value of P(N, /) is the coefficient of ¢/ in the expansion of 


77i(D;+ 4)! 


where q;=1—p,, and 7 runs over all units for which p; is not negligible. 
For three points on the lattice, namely on one of the lattice points, at the 
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centre of one of the triangular cells and at the centre of one of the sides of 
the cells, there is sufficient symmetry to allow computation of the 
coefficients up to f= 8, and taking into consideration the fifty units nearest 
to the axis. We find that the values of P(.V,/) vary significantly with the 
position of the axis only when P is very small. We have to bear in mind 
that in the integrals for which P(N,f) is required, the values of P corres- 
ponding to small values of V are weighted heavily. However, the values 
of P do not vary by more than 15% for f=3 and by smaller amounts for 
f >3, and we are justified in taking the mean values of P(N,f). This is 
confirmed by the measurements described later. Second, for / >8, we 
have considered only showers whose axes fall on a lattice point, and used a 
Monte Carlo method to determine the numbers of units triggered by 
showers of different V. From this data, we can determine P(N,/). 

From these values, we can also determine P(N, > f) the probability that 
a shower containing V particles should trigger at least f units. 

We assume for simplicity a spectrum of the form 


J(N)dN=kN~7" dN. 


Then the rate per unit area of showers triggering at least f units is 


7 00 


R(>f)= | J(N)AN. PW,> f) 


a” Riefy=k |" exp (—wy)P(w, >f)dw, where w=In VN. 
if we take a value yo, known nearly to satisfy the data and put y=y,)+6,, 
we may write 

R( >f=k [> exp (—wy)P(w, > f) de 


=e w exp (—wy)P(w, =f) ino+ 00? 
Both these integrals can be evaluated, and the different values of f from the 
experimental observations give simultaneous equations which can be 
solved by a least squares method for k and 5... We then calculate » 
from y=yo+6,. 

The spectrum calculated in this way refers to the number of showers 
striking the ground. It will be shown later (paper VI) that the solid angle 
for our apparatus is 0-5 radians. We can then calculate the spectrum 
per unit solid angle, and obtain 


Jl > N) = i LGN ON aa. 8-O 5 10 (105) NE eed? met hr-t sterc?. 
JN 
It is of interest to see how showers of different sizes contribute to the 
events which go into the various categories, i.e. the resolving power of 
the apparatus. In fig. 2 we have drawn dashed curves which show the 
resolution for the different categories. It will be seen that the resolution 
improves as the number of units struck increases, The resolution curves 
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are markedly asymmetrical, and in table 1 we give values of the mean of N, 
and N49, for which 10%, of the showers contain fewer particles than N49. 

The points marked on fig. 2 show the rates at which showers of different 
sizes are recorded. These are plotted at values of N given by 


9:5 x 10-9(108/NV):88 = | * I(N)P(N, =f) aN. 


Fig. 2 


Number of showers, m <hr 'ster | 


fe) 
1 


Table 1 
S38 =) >8 > 20 > 40 
N 5-5 x 108 V2 10? 2-710? 10108 | 2-8x 108 
Artes 3-9 106 | 9-5x106 | 18107 | 6-8x107 | 1-7> 108 
Nae 17108 | 473108 | 10x10? | 45%107 | 1-2 x 108 
Nir te 5:3 ><10® | 91110" 2-256 04 at-6 e 07 = 


This is the value of N at which an ideal apparatus, i.e. one with infinite 
resolving power would count showers at the same rate as the practical 
apparatus. The agreement of the ordinate with the value given by the 
line shows that the analysis is satisfactory, and that the assumption of a 
power law spectrum adequately accounts for the observations. 
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§ 4. THe MEASUREMENT OF NuMBERS OF PARTICLES 
The average number of particles in showers belonging to a particular 
size group has also been determined in the following way: , 
Let showers with V, N+dN particles be recorded in this group by the 
array at a rate f(N))dN per unit area. The average number of particles 
in this size group is then given by 


. | ‘ | ‘ 2m drp(r, N) f(N)dN 
| sean 


7cO 7c 
cage | Qnrdro(r,N)f(N)dN  . . . . . (1) 
PJoJo 
where P is the rate of detected showers in this size group and p(r, N) is the 
density of particles. 
The rate at which 7 detectors, each having an area S, will discharge in 
coincidence with showers in this group, is 


Gee if | FN) aN. 2nrdr{1—expl—p(r, NST} anc) 


Using the identity 


eqn. (1) may be written 


SNP= > . i | “ f(N) dN . 2ar dr[1 —exp{—pl(r, wysy¥} 
= 0 0 


L 


~s l co 1 , 
Therefore - xp( 55). reer neh wc ad 


For each size group N has been determined by recording R, from six 
G.M. counters, area 20cm?. The figures obtained are in the fourth row 
in table 1, where they can be compared with N obtained by using the 
equation 


ie J(N)P(N, > f)N aN 


| ” J(N)P(N, > f)dN 
0) 

It should be noted that in this measurement, apart from the requirement 
that the sum (3) should converge, no assumption has been made about the 
radial distribution of the showers. However, we have assumed that 
f(N) is independent of position on the array. This is equivalent to the 
assumption that the P(N,f) are independent of position as discussed 
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in §3. If f(N) varies with position, #, will be a function of position. 
To test this point, seven counters were placed at the points marked 1—7_ 
in fig. 3, inside a cell near the centre of the array. Thus R, which contri- 
butes most to the value of N was determined for each position. The 
values found are given in table 2, and are seen not to vary significantly. 

Later, the counters were all moved to the position marked 1 and the 
values for R, obtained are given in table 3. It can be seen that the 
convergence is satisfactory. 


Fig. 3 
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The errors in this method of determination can be estimated by the 
following considerations. The method consists in sampling the showers at 
different distances from the core. At large distances, the number of 
samples is large, and the density well determined. Moreover, the contri- 
bution to the total number of particles in the shower from large distances 
is small, so that the error introduced into the total is small. At short 
distances on the other hand, the number of samples is small, the density 
large, and not well determined by the six small counters, and the contri- 
bution to the total by particles from this region may be large. Thus for 
the >3 fold showers, the measurement of density becomes poor inside 
12 metres from the core. This region is sampled about 15 times during 
the experiment, and contains about 12% of the particles. The contri- 
bution to the error from this region is estimated to be about 8°%, and the 
total error on the number of particles should be about 10%. Similarly 
we find for > 8 fold showers an error of about 20°,. For > 20 fold showers, 
the region less than about 80 metres from the core triggers all the counters 
on nearly every occasion. ‘This region contains more than half the particles 


Observations on Extensive Air Showers: V 383 


in the shower. For these showers the method gives a low value for the 
number of particles, with a very large error. 

A further under-estimate of the number of particles in > 20 fold showers 
arises, since more than 5% of the showers falling near the boundary of the 
array discharge at least one of the small units but are not recorded as 
> 20 fold showers. 

Table 2 


Counter position 


Number of counts 


Ayal once dbittee Watch aaa pap, Dae 

27-9 18-0 14:8 13-1 11-5 123-5 

16-4 11-5 8-2 8-2 6:6 61-2 

9-8 6-6 4-9 3:3 3:3 30-9 

> 20 4-9 2-6 9-2 1-9 1-7 1-1 TD, 


§ 5. CONCLUSION 


We have measured the size spectrum of air showers containing between 
10° and 108 particles. We find 
is 


fC ZW y= 9-5 < 10-9(108/V)1:8920-12 m=" hr ster-*. 
This result may be compared with the M.I.T. experiment (Clark et al. 1957) 
who find 

JU 2 NN) = 8-2.x 10 4(109/ NV F806 m hr-*ster™. 
An independent measurement of the number of particles in the showers 
gives good agreement with these figures. 
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ABSTRACT 


Single crystals of nickel (purities 99-98% and 99-4%) have been deformed 

in tension at various temperatures between 4:2°K and 300°K. Shear 
stress-shear strain relations have been determined as well as the reversible 
change of flow stress with temperature. The work-hardening parameters 
obtained are discussed in terms of dislocation theory that has been developed 
mainly on the basis of data on Aland Cu. From the temperature dependence 
of the stress at the beginning of dynamical recovery the activation energy of 
cross-slip is estimated to be slightly lower than that of copper, implying a 
somewhat higher stacking fault energy. 
’ The advantages of the temperature change method for an investigation of 
the flow stress are pointed out following a theoretical analysis of the relation 
between this method and that of the work-hardening parameters. The 
present results on nickel confirm this relation. 

Metallographic and x-ray evidence has been obtained in favour of deforma- 
tion twinning in nickel at 4:2°k and 20°K although no drop in load has been 
observed of the kind found in copper. Serrations in the recorded load- 
extension curves indicate discontinuous flow at the lowest temperatures. 


§ 1. [INTRODUCTION 


AN understanding of the plasticity of metals with the help of dislocation 
theory seems to be possible so far only in relatively simple cases, such as 
single crystals of close-packed hexagonal and face-centred cubic metals. 
For the latter, the situation at low temperatures is particularly promising 
since, for example, three different mechanisms of work-hardening manifest 
themselves in separate stages of the stress-strain curve, as measured in 
a constant strain rate tensile test, of a single crystal (fig. 1). Following 
the Stuttgart school (Diehl! 1956 a, Seeger et al. 1957) the three stages can 
be characterized as follows: Stage I will be called the region of easy glide 
(Andrade and Henderson 1951). Here a small rate of work-hardening 
suggests that most of the dislocations produced during deformation leave 
the crystal through the surface. Stage LH, the rapid hardening range, 
is thought to be characterized by the formation of piled-up groups 
of dislocations behind Lomer—Cottrell dislocations. At high stresses, 
some of the slip dislocations become able to by-pass these obstacles 
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(‘eross-slip’). This process of dynamical recovery reduces the rate of 
work-hardening and leads to stage III of the stress-strain curve. 

Most of the earlier work on plastic deformation of face-centred cubic 
single crystals was done with aluminium at room temperature where, 
unfortunately, stage III sets in immediately (after a very short stage I). 
At 78°K, however, Al has been shown (Staubwasser 1954, Cottrell and 
Stokes 1955, Jaoul 1957) to exhibit stress-strain curves of the type 
drawn in fig. lf. This prototype has been established mainly by tensile 
tests of copper crystals at 300°K and below (Diehl 1956 a, Blewitt 1953, 
Blewitt et al. 1955, Adams and Cottrell 1955, Rebstock 1957, Berner 1957, 
Andrade and Aboay 1957). Other experimental evidence, especially 
electron-microscopic observations of slip lines (Diehl et al. 1955, Mader 


Fig. 1 


shear stress 


Pl 
shear strain 


Schematic shear stress-shear strain curve of face-centred cubic single crystal at 
low temperatures. 


1957), has helped to connect certain dislocation processes with the various 
stages of the stress-strain curve, as indicated above. There is, however, 
a great need for experimental data on work-hardening ‘parameters ’ 
—characterizing quantitatively each of the three stages—obtained from 
a single set of crystals, at different temperatures. It was thought that 
a third face-centred cubic metal such as nickel would offer the necessary 
check on theories developed on the basis of work on Cu and Al. As will 
be seen below, the extent of the three stages changes strikingly with 
temperature between 300°K and 20°K, even for a high melting point 
metal such as nickel. 

The measured values of the parameters are given in §3, and in §5 some 
of the measurements, especially on dynamical recovery, are evaluated 
in terms of existing dislocation theory. In addition to the temperature 


+ There are also some published data on the first two stages of work-hardening 
of Al single crystals at 4-2°K (Sosin and Koehler 1956, Noggle and Koehler 
1957). 
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variation it is now well known from work on Al (e.g. Staubwasser 1954, 
Lange and Liicke 1952, 1953, Davis ef al. 1957) and Cu (Diehl 1956 a, 
Rosi 1954, Paterson 1955) that the shear stress-shear strain curve still 
depends on the orientation of the crystal with respect to the tensile 
direction. For this reason crystals of one orientation were investigated 
at different temperatures as well as different orientations at a given 
temperature. 

The influence of crystal orientation makes it difficult to evaluate previous 
data by Andrade and Henderson (1951) on work-hardening of Ni. These 
authors measured stress-strain curves of four crystals at 90°, 300°, 508° 
and 833°K by a stepwise loading technique (99-99, Ni, unknown 
orientation). 

In addition to the type of analysis of work-hardening described above, 
Cottrell and co-workers (Cottrell and Stokes 1955, Adams and Cottrell 
1955) have introduced a method of separating the flow stress of a cold- 
worked crystal into two contributions, one due to the dislocation structure 
and another due to the dislocation mobility. They measured the reversible 
change in flow stress Av connected with a change in temperature (7', — 7’3) 
at a given strain. Az determines the change in dislocation mobility, or 
dislocation friction, with temperature in a given dislocation structure. 
We will give our results on A7(7',) for nickel in§4. In $5 we will discuss 
the relation between the two methods of analysis of work-hardening and 
their results. 

In recent papers Blewitt (1953), Blewitt et al. (1955, 1957) established 
evidence for discontinuous slip and deformation twinning in copper 
single crystals at 4°K and 78°K under very high stresses. We have searched 
for such a phenomenon in Ni and Cu at 20°K (and in one example, using 
Blewitt’s equipment, also at 4°K). The results of these tests are given 
in $6. 

§ 2. EXPERIMENTAL PROCEDURES 
2.1. Specimen Preparation 

Single crystals of nickel were grown from the melt under high vacuum 
in pure Alumina crucibles by the Bridgman technique. Most of the 
crystals were of Johnson Matthey J.M. 890 nickel, nominal purity 99-999%. 
Spectroscopic analysis at the Chicago Institute before and after crystal 
fabrication indicated a purity greater than 99-98%. <A gas analysis by 
vacuum fusion after different heat treatments of the crystals showed 
about 10% N,, 5x 10-°% H, and 10-°% O,. A micro-carbon analysis 
gave less than 10-°% C for vacuum annealed crystals (about 3 x 10-2% C 
for the hydrogen pre-annealed crystals). Some crystals, marked ‘c¢’ in 
front of the number, were of commercial purity, 99-4°% Ni, with Co and Fe 
as Major impurities. 

The dimensions of the as-grown single crystal wires were as follows: 
diameter 0-088 in. and maximum length 8 in. After growth most of the 
crystals adhered to the crucible walls by a thin blue layer of what is believed 
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to be a mixed Ni-Al oxide. The best way of removing the crucibles was 
with the help of a micro-sandblasting device, using Al,O, powder. The 
thin jet cuts the brittle crucible into small pieces which fall off without 
damaging the soft nickel crystal. Following this the crystals were cleaned 
in HNO, and were cut into 2-8 in. long pieces by means of an acid saw 
as described by Maddin (1950). In this way a number of equi-oriented 
specimens grown under the same conditions could be investigated. 

After preparation the crystals were usually annealed in vacuo for 12 hours 
at 900°C, The mechanical behaviour of some as-grown crystals and of 
some crystals that were annealed for 12 hours in wet or dry hydrogen 
before the vacuum anneal was compared with that of vacuum annealed 
crystals of the same orientation. As will be seen in §3, the as-grown 
crystals showed a higher critical shear stress 7), and the hydrogen 
pre-annealed ones a lower 7) than the standard vacuum annealed crystals. 


Fig. 2 


Stereographic plot of the orientation of the tensile axes of nickel single crystals. 


After annealing, the x-ray photographs showed sharp Laue spots. 
The orientations of the wire axes were determined by a simple rotating 
crystal method (Schmid and Boas 1935); they are plotted in the 
stereographic unit triangle in fig. 2. The crystals were subsequently 
soldered, using a micro-flame, into steel grips, the ball-shaped ends of 
which fitted into the holders (IL in fig. 3) of the tensile machine. At first 
crystals to be tested at low stresses were soft soldered. Later a special 
low melting silver solder was used with little effect on the critical shear 
stress. A more complex holding device employing a piece of piano wire 
for better initial alignment was also tried and did not have much effect 
on the critical shear stress. Finally, the gripped crystals were chemically 
polished. 

2.2. Method of Deformation 

The crystals were strained in a Tinius Olsen tensile testing machine 
(employing the ranges 100, 150, 200, 1000 lb.) which had been adapted 
for use at temperatures down to 20°K. The low temperature unit Is drawn 
in fig. 3. In this arrangement the crystal is pulled between the stainless 
steel rod J and the base N which is kept at a fixed distance from the top 
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Fig. 3 
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Low temperature extensometer unit (described in the text). 
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plate B by four stainless steel rods M. The pull rod J is connected to the 
pendulum system of the machine, and the top plate B to the moving 
_ cross-head, by means of the clamps A. A Dewar flask F surrounds the 
assembly. It is fixed to the top plate B by a brass case E having slot 
windows. The case can be made gas-tight for the tests under liquid 
_ hydrogen by means of Teflon O-ring seals O, C and K and by connecting 
it to an exhaust line at D. A piece of styrofoam G and an outer N, gas jet 
prevent the moving seal K from freezing after liquid hydrogen is transferred 
into the Dewar through C. A filling of 6 litres of liquid hydrogen was 
sufficient for an hour’s run. Further tests were made using liquid nitrogen 
(78°K), a solid CO,-trichlorethylene mixture (~200°K), and at room 
temperature. 

The gauge length of the tensile specimen measured about 2in., and the 
cross-head speed was kept constant at 0-05 in./min giving a strain rate 
between 10~* and 10~°/sec depending somewhat on the load range 
because of the softness of the machine. Strain was measured by a linear 
transformer between B and J, fig. 3, and was continually recorded against 
load. The strain measurement has been checked by microscopic observation 
of the distance of marks on the crystal before and after deformation. 


§ 3. STRESS-STRAIN CURVES AND WoRK-HARDENING PARAMETERS 


From the recorded load—extension curves shear stress-shear strain 
curves were computed by using the well-known relations (Schmid and 
Boas 1935). The initial orientation determines the position of the tensile 
axis in the dodecahedral plane through [111] and [100] in the stereographic 
triangle where double slip sets in. (Not much overshooting has been 
observed (Osswald 1933).) From thereon the vonGoler-Sachs (1927) 
method for double glide has been used to compute shear stress and shear 
strain. It gives a considerably higher work-hardening rate than the 
application of the single slip formula which is here unjustifiedt. Typical 
stress-strain curves of nickel at different temperatures and for different 
orientations are given in figs. 4 to 10. The appearance of some of the 
~ curves is somewhat complicated by the temperature changes introduced 
to determine the reversible change of flow stress. Fortunately, after a 
few percent strain past the point of the temperature change the work- 
hardening parameters settled down to the values characteristic of the new 
temperature—except during stage I (discussed below). The parameters 
defined in fig. 11 have been obtained from the stress-strain curves. The 
results will be discussed together below in terms of certain common 
parameters evaluated from ail crystals. In this way the variation with 
temperature and orientation can be best recognized and separated from 
the individual scatter between the crystalst. 


+ Apparently, the curves given by Andrade and Henderson (1951) and 
Andrade and Aboav (1957) have been calculated for single slip throughout. 

+ Also, due to the temperature changes, not all the parameters could be taken 
from every crystal, 
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Shear stress-shear strain curves of high purity nickel crystal 6. 
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Shear stress-shear strain curves of high purity nickel crystal 14. 
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Shear stress-shear strain curves of high purity nickel crystals 22c and 31. 
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Shear stress-shear strain curves of commercial nickel crystal ¢19. 
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Definition of the work-hardening parameters. 


3.1. The Critical Shear Stress ro 


This is defined by the intersection of the extrapolated linear easy glide 
region with the stress axis. Table 1 gives the measured values. An 
increase of rt, at lower temperatures is evident, but no orientation 
dependence (see also Diehl 1956a), within the considerable scatter of 
this quantity. A better way to determine the temperature variation 
of 7, for one crystal is by studying temperature changes during stage I, 
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after a few per cent strain. The results described in §4 indicate a value 
of Ar,=150g/mm2 between 78°K and 300°K, and A7ty=350g/mm? 
between 20°K and 300°K, This is consistent with the average values 
from table L: 


7) = 600-850 g/mm? at 300°K. 
Ty) = 750-850 g/mm? at 78°K, 
T) = 900-1100 g/mm? at 20°K. 


Special attention is called to the hydrogen pre-annealed crystals 5a and 7b 
which have critical shear stresses about 230 g/mm/? lower than the standard 
vacuum annealed crystals. As-grown crystals (5b, 7¢ and ¢19a) show 
up to 50% higher +, values than the vacuum annealed crystals. No 
particular difference is visible between 99-4 and 99-98 nickel regarding 79 
at 300°K. 

Table 1. 7, (kg/mm?) 


Parameter Values at Temperature 
Crystal ; 
: Remarks 


No. ‘ 78°K 200°K 300°K 


0-80 0-65 
*mishandled *(1-05) 0-85 0-85 
0-75 0-65 


0-60 

0-85 D8: 0-75 

0-57 

pre-annealed in H, 0-33 

as grown 0-75 
2 


pre-annealed in H, 0-47 
as grown 1-05 


0-70 


4-2°K : 1-15 0-80 
0-60 
as grown 1-15 
0-75 


Andrade and Henderson (1951) give t)=1040g/mm? at 300°K, 
Ty = 1360 g/mm?* at 93°K for their 99-9 nickel. Osswald (1933) measured 
an average 7)=580g/mm? on four 99-79% nickel crystals at room 
temperature. 


3.2. The Work-Hardening Rate },; during Easy Glide 


This has been compiled in table 2. There seems to be a slight decrease 
of 3; on changing the temperature from 300°K to 78°K and 20°K within 
the uncertainty of the scatter to which 9, is particularly susceptible 
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due to growth conditions (Diehl 1956 b). 9; increases strongly (from 3 
to 12kg/mm®) for unfavourable orientations where other slip systems 
come into play early. As will be shown below, the extent of stage I in 
such cases is very small so that 9, might be influenced by the transition 
into stage II. The relatively small 9, of the hydrogen pre-annealed and 
commercially pure nickel crystals, especially the unfavourably oriented 
crystal ¢ 19, should be noted. We do not know the reason for the abnormal 
easy glide parameters of the crystals pre-annealed in hydrogen instead of 
vacuum (see also Cupp and Chalmers 1954). Hydrogen was shown not 
to be responsible for the yield phenomenon observed (Haasen and Kelly 
1957). 
Table 2. 9, (kg/mm?) 


Parameter Values at Temperature 


Sine Remarks 
; 20°K 78°K 200°K 300°K 
6 3°25 2-75 4-0 
l4 3-0 7-5 
31 3:6 5:5 4:5 
22a 6-2 
22c 10-0 
18 5-0 75 
Ta 12-0 
7b | pre-annealed in H, (aS) 
7¢@ as grown 8-5 
11 11-5 11-0 
5a | pre-annealed in H, 7:3 
5b | as grown 
oc 
16 
12 4-2°K : 12-5 
e23 8:5 
e19b 5 


3.3. The Extent of Hasy Glide ay, (Table 3) 


This is strongly temperature dependent as already noted by Andrade 
and Henderson (1951). In our crystals a,;; increases from 7% at 300°K to 
42°, at 20°K for a favourable orientation. For unfavourable orientations 
a;; is much smaller and varies less with temperature. Commercial nickel 
shows a different behaviour. Here a,, is relatively large even for the 
unfavourably oriented crystal ¢ 19 (fig. 10) and seems to be independent 
of temperature. . 

Figure 12 summarizes a,,(7') values for pure and commercial nickel. 
When easy glide begins at a higher temperature 7’, and continues at a 
lower temperature 7’, the large a,;; associated with 7 is obtained only if 
the glide at 7’, had not yet touched the transition into stage Il (compare 
erystal 12a (fig. 9), 14¢ (fig. 5), 22 ¢ (fig. 6) with 31 b (fig. 6), ¢ 9a (fig. 10)). 
Thus once rapid hardening had started at one temperature the crystal was 


2E2 
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never able to continue easy glide at any other temperature. This behaviour 
differs from that at the transition from stage II into stage III where 
crystals pulled into stage IIL at a higher temperature came back into 


Table 3. ayy 


Crystal Parameter Values at Temperature 
No. Remarks 


20°K 78°K 200°K 300°K 


6 0-42 0-30 0-075 


14 0-17 
31 0-35 | 0-13 (0-10) 


0-15 


0-095 


pre-annealed in H, 
as grown 


pre-annealed in H, 


5b | as grown (0-025) 
5e (0-05) 
16 (0-05) 
Ly 4-2°K : 0-11 
c23 0-08 
b 0-10 


as grown. 


Fig. 12 
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Extent of easy glide a;,; vs temperature for nickel crystals of different purity 
oriented favourably (max) or unfavourably (min) for single slip. 
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stage IT at a lower temperature. This difference is in agreement with the 
dislocation mechanisms believed to be operating in stages I and II 
(see § 5). 


3.4. The Stress at the Setting-in of Rapid Hardening, r,, (Table 4) 

This follows from 9, ty and a;; according to 9, Y(r;;—T)/ayy. The 
increase of 7; by a factor of 2-5 connected with a change in temperature 
from 300°K to 20°K for a favourable} orientation reflects the increase of a, 
under the same conditions since 9, and z, are less temperature dependent. 
Tz; Seems to be smaller for unfavourable orientations} at 20°K and 78°K. 
At 300°K, however, 7,; is somewhat larger for the unfavourable 
orientations}. Therefore the effect of temperature on 7,, is reduced for 
these orientations as it is on ay. If one considers 7,;; and 9, as primary 
quantities responsible for the extent a, of easy glide then the temperature 
dependence of aj; is largely due to that of +,,; while the orientation 
dependence of a;; is mainly that of 9. 


Table 4. 7,; (kg/mm?) 


Pec ctal Parameter Values at Temperature 
No Remarks 
; 20°K 78°K 200°K 300°K 
) 6 2-5 1-6 0-95 
14 1-55 
31 2:3 1-3 
22a 1-8 
22¢ 0:95 
18 1-3 1-25 
Ta 1-15 
7b | pre-annealed in H, 090 | 
7c as grown 125 | 
ll a | 
5a | pre-annealed in H, 0:80 | 
5b | as grown (<1:35) | 
5e (<1:5) 
16 1-95 | 
12 4:2°K : 2-5 
c23 1-3 
cl9 (2-9) 1-25 
as grown 


+ We continue to call unfavourable orientations those near the [111]-[100] 
boundary of the stereographic triangle, favourable those far from it which 
ullow considerable single slip to occur in the crystal. 

tIt has been proposed by Rosi (1954) and Garstone, Honeycombe and 
Freetham (1956) that (7,;-79) OF 71/7 is a more or less orientation-independent 
constant quantity determining the beginning of stage II similarly to Schmid’s 
aw for ty. This does not hold for our crystals (neither did it in Diehl’s investi- 
sation of Cu). Also, r;; varies more strongly with temperature than does Tp. 
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3.5. The Slope of the Rapid Hardening Range, >, 


This is shown in table 5 and fig. 13. It shows a stronger temperature 
dependence than expected from the measurements on copper made by 
Blewitt (1953), Blewitt et al. (1955), Berner (1957) and Andrade and 
Aboav (1957). The dependence on orientation lies within the scatter 
between individual crystals. No difference between 99-98 and 99-4 nicke 
has been found. Taking weighted averages from table 5 and fig. 13 we 
consider as best values for 9;,;: at 20°K, 28kg/mm?; at 78°K, 25 kg/mm?; 
at 300°K, 21kg/mm2. Osswald (1933) observed for a 99-7% Ni at room 
temperature a value 9, = 26 kg/mm’. 


3.6. The Stress at the Beginning of Dynamical Recovery, 74, 


This is a very informative parameter of work-hardening. The measuredt 
values in table 6 show that 7,,; increases strongly at lower temperatures. 
No orientation dependence can be recognized. For the commercial nickel 
crystal ¢ 23 7,;; is 20% larger than for the pure nickel crystal 18 of the ~ 
same orientation. As will be discussed below, theoretical calculations — 
predict a relation of the form log7,;,;;=A— BT between the stress at the ; 
beginning of stage III and the temperature 7’ of deformation. A logarithmic 
plot of 7,;,; against temperature does indeed give a straight line (fig. 14) — 
with the constants being B=4-4 x 10-3, A=2-98 when 7,,;; is measured — 
in kg/mm? and 7’ in °K. We have also plotted 7,;, for Cu from the data 
of Blewitt et al. (1955) (crystals nos. 304 and 305) and Berner (1957). 
Again the logarithmic relation holds well, the constants being B= 4-1 x 10-3 
and A=2-48 (Blewitt), B=4-4x10-3 and A=2-72 (Berner). We will 
use these results in §5 to estimate the activation energy of cross-slip in Ni 
relative to that of Cu. 


3.7. The Work-Hardening Rate during Stage I11, 3; (Table 7) 


This depends on temperature in much the same manner as %y,, the 
ratio Dy;_/9;, being constant, about 0-60 for 78°K, 200°K and 300°K; there 
is not much dependence upon orientation. However, 9,,;/9;; increases 
to about 0-75 at 20°K (and 4°K) with considerable scatter from crystal to 
crystal. The linear character of stage III is by no means as pronounced 
as that of stage I and stage Il. The work-hardening rate drops below 9444 
at the end of the stress-strain curve so that $,,;; has in some cases the 
meaning of a tangent which intersects the stage II line at 7,,,; as defined 
above. 


} Since stage ITT is not always linear 7,,; has been defined as the stress on the 
stage I line where the measured curve departs by 200 g/mm?. 7; so defined 
agrees with the stress at the intersection of the stage Hf and stage UT straight 
lines where the latter is well developed. 
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Fig. 13 
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Rate 3,, of rapid hardening vs temperature. The crosses mark measured best 


values, the circle at 20°K the value expected from temperature change 
experiments (§ 5). 


Fig. 14 
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Stress 7;;; at the beginning of dynamical recovery, measured (in kg/mm?) on 
nickel and copper (Blewitt e¢ al. 1955, Berner 1957), plotted logarith- 
mically vs temperature. (The crosses mark averages from table 6 
including the high values.) 
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Table 5. 9%, (kg/mm?) 


Parameter Values at Temperature 


Crystal emurke ge at ke sgh ee 
Bes 20°K 78x | 200% | 300% 
6 30-0 27:5 20-5 

14 27-5 22:5 
31 30-0 22-0 23-0 
22a (28-0) 
22¢ 28-5 22-5 
18 23-0 20:3 20-0 19-5 
Ta 22:0 22-0 
7b | pre-annealed in H, (24-0) 
7c | aS grown 24-0 20-0 
11 26-0 27-5 
5a | pre-annealed in H, 26-5 
5b | as grown 29-5 21-5 
5e 26-0 23-5 
16 27-0 20-0 
12 4-2°K: 31-5 32-0 25-5 22-0 


Table 6. ty, (kg/mm?) 
Grell . | Parameter Values at Temperature ; 
No emarks ; 
; 20°K 78°K 200°K 300°K 
6 17-0 15-2 6-2 
14 Ie 14-7 
31 17:0 (16-0) 4-9 
22a Lie 
22¢ 20-3 
18 15-4 8:7 7:6 
7a ¢ 
7b | pre-annealed in H, q 
7¢ as grown Q 
11 (13-3) 
5a | pre-annealed in H, 20-7 
5b | as grown 21-0 12-7 
5c (22-0) 5-5 
16 2 
12 4-2°K : 21-2 18-1 
e23 18-5 
cl9 2 
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Table 7. dy; (kg/mm?) 


Parameter Values at Temperature 
Remarks 


20°K 78°K 200°K 300°K. 


pre-annealed in H, g 
as grown 2 


pre-annealed in H, 23-5 


4-2°K : 23-0 


Table 8.  ayyz (@,) 


Crystal Parameter Values at Temperature 


No Remarks ooo 
‘ 20°K 78°K 200°K 300°K 
6 0-92 (1-32) | 0-80 (1-10) 0-33 (0-86) 
14 0-80 (1-28) | 0-75 (1-12) 
3l —(1-20) 
0-84 (1-15) | 0-72 (1-08) 0-30 (—) 
22a 0-68 (0-97) 
22¢ 0-72 (0-98) 
18 0-70 (0-96) | 0-42 (0-88) | — (0-80) 
0-35 (—) 
7a 
7b 
7¢c 
11 0-54 (0-90) 
5a | pre-annealed in H, | 0-74 (0-98) 
5b | as grown 0-67 (0-85) | 0-40 (—) 
5¢ 0-81 (1-12) 0-22 (—) 
16 
12 4-2°K : 0-71 (1:04) | 0-57 (0-88) 
c23 0-80 (0-98) 
cl9 
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3.8. The Strains at the Beginning and the End of Stage III, 
Ay; and a, Respectively 


These are given in table 8. They depend considerably on the deforma- 
tion-temperature history of the crystal. Of particular importance for 
these quantities is the extent of easy glide a;, as determined by the 
orientation of the crystal and the temperature during the stage I 
deformation. We have seen ($3.3) that a ,; increases strongly at low 
temperatures, and for favourable orientations. The orientation dependence 
of ay;, and a, is essentially that of a; (see the preceding sections) and is 
therefore determined by the variation of 9; with orientation. Concerning 
the temperature influence on a,,; and a, a contribution besides that of aj; 
comes from 7;;;(7'). The extent of stage II as measured by (a;,;;— yy) 18 
about 50 to 60%, for low temperatures and about 20° at 300°K for this 
reason. Altogether, a,,; increases by a factor of three at low temperatures 
for a favourable orientation. The extent of stage III, (a,—a,;,), is about 
40% +10% at all temperatures investigated. Therefore the strain at 
fracture, a,, increases from about 85°% at 300°K to 130% at 20°K. 


3.9. The Stress at Fracture, 7, 


This is given in table 9. It is independent of orientation within the 
scatter. Its temperature dependence is predominantly that of 7,,,; although 
(7,—T 1) increases somewhat between 78°K and 20°K. The stress at 
fracture at 20°K is about G/400, where G is the shear modulus. 


Table 9. 7, (kg/mm?) 


Crystal Pe : Parameter Values at Temperature 
NO emarks 
20°K 78°K 200°K 300°K 
6 24-6 19-5 12-2 
14. 31-8 20-4 
31 23°7 23-0 20:8 
22a 23-0 
22¢ 26-0 
18 18-7 13-5 11-1 
7a 
7b 
7¢ 
iil 18-8 
5a | pre-annealed in H, 25:5 
5b | as grown 24-5 
5¢ 28-0 
16 
12 4-2°K ; 27-8 24.3 
e23 21-0 
cl9 
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§ 4. THe REVERSIBLE CHANGE Or FLow Srress wrrn TEMPERATURE 

The advantage of employing the method of temperature changes during 
deformation has already been indicated above. It enables one to determine 
the change of flow stress with temperature in a given cold-worked state. 
The experimental procedure is that one deforms a crystal at a temperature 
T, up to a stress 7(7’,), unloads, changes to 7’, and resumes deformation 
at 7’, which starts now at a stress 7(7',).  (7(7',)—7(7',)) can thus be 
determined as a function of the temperature 7’, and the work-hardening 
(measured by 7(7',)t), 7, being fixed at a low temperature. Typical 
examples of the change in the flow stress (7, —7,) of nickel with temperature 
are shown in the stress-strain curves of figs. 5, 6, 7, 8, 9, 10. In order to 
rationalize the results of such experiments Cottrell and Stokes (1955) 
and Adams and Cottrell (1955) in previous investigations on Al and Cu 
have plotted the flow stress ratio 7(7',)/7(7,) = 1 — (7(7,) —7(1g))/7(7,) as 
a function of 7’, for temperature changes made at large strains (7, < 7’,). 
In this way they take advantage of two observations: (i) that 7/7, 
becomes independent of the strain a for a> 10-15%, which is in stage III 
at the temperatures 7’, they used, and (ii) that the stress ratios can be 
determined by stepwise temperature changes according to the relation 
T/T, = (79/T3)(73/7,). Since they found that the cold-worked state produced 
at the low temperature was unstable at high temperatures for large strains 
( work-softening ’) Cottrell and his co-workers used only the high-to-low 
temperature transitions. 

In the present work it was found useful (see also Seeger et al. 1957) to 
evaluate the data by a plot of (7,—7,) against 7, which gives the linear 
relation 

Ty —T2= Ayo t Oy a7 
already at relatively small strains, i.e. in stage II (not, however, during 
easy glide). This relation turns into Cottrell’s 7,/7,=constant for large 
stresses T, >, where 7,/7,=1/(b;.+1). Figures 15 to 18 show the results 
for (7,;—7T.) as a function of 7, on nickel for the temperatures 7’, and 7’, 
noted{. The most detailed results, those on the temperature changes 
between 78° and 300°K (fig. 15), give an idea that the linear relation holds 
within reasonable limits for different crystals and different orientations. 
It appears that in this way the scatter is averaged out better than by 
determining individual ratios of 7,/7, at large 7,. One can see that 
commercially pure nickel crystals show an anomalously large change in flow 
stress with temperature during easy glide (figs. 7 and 10). Even for pure 
nickel the (+, — 7.) values determined during easy glide lie consistently above 
the line obtained from stages Il and III. The constant value of (7, —75) for 


‘small 7, fits with the one determined from the temperature dependence 
OS ee eee re 

+ 7, denotes 7(7',) in this section. 

+ We have always subtracted the transient yield phenomenon. according to 
Haasen and Kelly (1957) in the determination of (7-74), le. we have used the 
lower yield stress at the new temperature. This was not consistently done so by 
Adams and Cottrell (1955) and Rebstock (1957). 
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of the critical shear stress ($3). We should point out that temperature 
changes in both directions lead to about the same change in flow stress, 


and that no actual drop in load has been observed at 300°K after 
deformation at a lower temperature, as occurred in Cu and Al. 


Fig. 15 


res i ah a 


Change in flow stress after a temperature change between 78°K and 300°K, 
measured at different values of the room temperature flow stress T5399. 
(The arrows give the direction of the temperature change.) 
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Change in flow stress after a temperature change between 20°K and 300°K 
measured at different values of the room temperature flow 
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From figs. 15 to 18 we get for the changes between the temperatures 
of liquid hydrogen (H), nitrogen (N), dry ice in trichlorethylene (1) and 
room temperature (R) the values by), = 0-27, byy=0-065, Dyy=0-167 or 


Fig. 17 


Tg" Tq (kg/mm?) 


Te(kg/mm?) 


Change in flow stress after a temperature change between 20°K and 78°K 
measured at different values of the flow stress at 78°K, T7.. 


Fig. 18 


T309 (kg/mm?) 


Change in flow stress after a temperature change between 200°K and 300°K 
measured at different values of the room temperature flow stress T3099. 


0-18 and bp;=0-05. From these 6’s the ratios 7)/rty=1/(by_ +1) were 
calculated. For the transition between 20° and 300°K the flow stress 
ratio could also be obtained from 


T390/T20 = (390/778)(T78/T29) = 97805 or 0:796 


compared with 7399/T2)= 0-788 by direct transition. 


Ta90/ UP ee (T290/ 7390) (T300/ To) = 0-827 
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was obtained only via the transition to room temperature. The flow stress 
at 7’, relative to that at 7’, = 20°K has been plotted against temperature 7’, 
in fig. 19. For comparison 7»/t,) for Cu and Al are also given in fig. 19. 
Values of 7/75) for Cu with 7 between 473°K and 90°K, measured by 
Adams and Cottrell (1955), were referred to 20°K according to 
T/T a9 = (T7/Te9)(T99/T20) With Ty9/7299 = 0-978 taken from our own measure- 
ments on Cu crystals. For Al the data of Cottrell and Stokes (1955) 
between 78°K and 423°K were combined with those of Basinski (1957) 
on polyerystalline, solution-treated 248 Al between 20°K and 80°K. They 
fit together surprisingly well. Comparing the three face-centred cubic 
metals in fig. 19, the curve of nickel is intermediate between those of Cu 
and Al and on the whole more similar in shape to that of Al although 
the change in slope occurs at a somewhat higher temperature (190°K for 
Ni, 120°K for Al). 


Fig. 19 
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Ratio of the flow stresses at temperature 7' and at 20°K vs temperature 7’ for 
nickel, copper (Adams and Cottrell 1955), and aluminium (Cottrell and 
Stokes 1955, Basinski 1957). (The average values of Adams and 
Cottrell’s measurements on Cu have been referred to 20°K as stated 
in the text. — T399/Tg9 is taken from our own measurements on copper.) 


§ 5. Discusston 


It is now generally agreed that work-hardening in face-centred cubic 
crystals is caused by two types of dislocation processes. In the first 
process, the dislocations spreading slip along the main glide planes cut 
other dislocations which do not lie in this plane. The production of jogs 
in this process and the movement of such jogs in screw dislocations by the 
creation of point defects require extra work to be done by the applied 
stress. In the second process dislocations are stopped in the crystal by 
reactions with dislocations of other slip systems leading to the formation 
of Lomer—Cottrell dislocations. More dislocations are piled up behind these 


— ee 
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obstacles, and the elastic stress fields arising from these groups work- 
harden the crystal. In the following discussion these two Sonu nnions 
to the flow stress will be called (cf. Seeger et al. 1957) ry (‘ Schneidprozesse , 
dislocation cutting) and rt, (the elastic interaction of dislocations ie 
proportional to a shear modulus G’). We will now discuss our measure- 
ments on work-hardening and flow stress during stages I and IIL in 
terms of these dislocation processes (§5.1). Stage III, the ‘ dynamical 
recovery’ region (Diehl et al. 1955, Seeger et al. 1957), will be dealt with 
separately in $5.2. We will base this discussion on recent reviews of 
the present state of the dislocation theory of work-hardening (Seeger 
et al. 1957, Seeger 1957) which was developed mainly on the basis of 
experiments on Al and Cu. The measurements on Ni support and 
extend these conclusions in many respects. In §5.1 we will stress in 
particular the relation between the two methods of analysis of work- 
hardening used (§§3 and 4) which became evident during this work. 
The result of this discussion makes the temperature change method 
appear even more useful and attractive. In $5.2 an evaluation of the 
data on the setting-in stress of stage III as a function of temperature 
leads to an estimate of the stacking fault energy of Ni relative to that 
of Cu. 


5.1. Stages I and IT 


During stage I the crystal slides essentially on one set of slip planes. 
This means that only few Lomer—Cottrell dislocations are formed and 
the r,, contribution is small. The long slip lines (Diehl et al. 1955, Mader 
1957) and the dependence of stage I on the specimen diameter (in the 
order of magnitude of 1 mm) (Rebstock 1957) show that many dislocations 
leave the crystal through the surface. Accordingly, the rate of work- 
hardening is low as it is in close-packed hexagonal crystals. This ‘laminar’ 
flow is very sensitive to disturbances introduced by inhomogeneity of the 
deformation, growth defects of the crystal, and also the interference of 
other slip systems which is likely to occur at this stage if the orientation 
is unfavourable (near the boundaries of the stereographic triangle (Haasen 
and Leibfried 1952, Haasen 1953)). Accordingly, the work-hardening 
rate, 9,, during easy glide depends strongly on orientation in our measure- 
ments on nickel crystals as well as those on copper. 

The temperature dependence of the flow stress in this first stage is 
about that of the critical shear stress which has been treated quantitatively 
by Seeger (1955): the contribution of the dislocation cutting processes 
depends on temperature in the form 


U —kT log (é/€) re ae 
EY fp Se NE RL T<T,= ———-~ 
rst) ; 2 ° Flog (Je) ; (1) 
| r.(Z') = 0 for T2T, | 


U 
where U is the activation energy of the process (jog formation or creation 


of vacancies and interstitials), v=bdl, b=strength of the dislocation, 
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—— ee 


d=diameter of the intersected dislocation or of a point defect, /= spacing — 


of the dislocation forest or distance between jogs, é=strain rate. The 


flow stress contribution tg of the elastic interaction of dislocations — 


in the crystal depends on temperature only as the shear modulus: 
tTg=abG4/N where N is the dislocation density (cm~*), « a constant of 
order of magnitude 1/10 (Seeger 1957). 

One thus obtains a critical shear stress t)=7g+7q Which decreases 
linearly with increasing temperature at low temperatures and is about 
constant at higher temperatures. Our results on 7,(7') by direct measure- 
ment on different crystals or—more reliably—by temperature change 
during stage I agree with such a general behaviour, the main rise in 7» 
for nickel being below 100°K. A quantitative analysis of 7, in terms of 
eqn. (1), as in the case of copper (Seeger 1955), faces serious difficulties : 
in order to separate 7, one has to know the temperature dependence of 
the shear modulus of as-grown nickel crystals. This depends sensitively 
on crystal perfection, due to the ferromagnetic AH-effect, even below 
room temperature. It is believed that the large scatter we observed in 
the critical shear stress just reflects different elastic moduli. After some 
percent deformation the AH-effect almost disappears and a unique 
temperature dependence of the moduli may be determined (see p. 410). 
Thus the method of temperature changes after some deformation, as 
described below, seems to be more reliable for an analysis of 7,(7') than 
the direct measurement of the critical shear stress, in the case of nickel 
crystals. It should be noted, though, that the measurements of Andrade 
and Henderson (1951) not only show higher absolute values of +, but 
also a 30° increase in 7, going from 300°K to 90°K while ours exhibit 
only a 10% change in the same temperature range. 

The stress 7;; where rapid hardening starts is thought to be identical 
with the critical shear stress of certain secondary glide systems (Seeger 
1957), the dislocations of which efficiently block the path of dislocations 
of the main slip system. From eqn. (1) it follows that the temperature 
dependence of 7,; is stronger than that of t, since v~, i.e. the density of 
the dislocation forest, is larger as seen from the secondary systems (many 
dislocations of the main system cutting) than seen from the main system 
(fewer dislocations of the secondary systems intersecting), see also Barrett 
(1952). This strong temperature dependence of 7,; explains the strong 
rise of a; at lower temperatures. 

After the stress r,; is reached the secondary slip systems continually 
interfere (Rebstock 1957) with slip on the main system which carries the 
deformation according to the law of maximum shear stress. No break 
in the shear stress-shear strain curve could be observed when the tensile 
axis reached the symmetric position where the conjugate slip system 
takes over (Osswald 1933). During stage IT the length of the active slip 
lines has been observed (Rebstock 1957) to decrease with increasing strain 
which is consistent with the idea of a formation of more and more 
Lomer—Cottrell dislocations. The groups of dislocations piled up at these 
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obstacles make the 7, contribution to work-hardening gain in importance. 
ty is, however, also rising with deformation as v-! increases. It will be 
‘shown now that the temperature dependence of 9,,, the rate of rapid 
hardening, and of 7,/7,, the reversible change in flow stress, is largely due 
to 7,(7’). Furthermore, both these quantities give identically the same 
‘information on 7,. 
_ We start from our empirical result that during stage II (and III) the 
flow stress before and after a change in temperature (between the 
temperatures 7’, and 7',) depends on work-hardening (as measured by 75) 
in the form 

Ty —T2=4+ dry. NE eS AND 5 Tae Anny SE) 


The flow stress at any temperature is composed of the two contributions 


T=TytT¢ Sp pi ete see) 
where for a given dislocated state of the crystal the definition of 7, implies 
TE Gy 

==. west Aeey fA ed en: 
Ta, Ge or 
Accordingly (2) turns into 
G G 
ee +n(Z=1), =e eyo s gh?) 


Now 7, also depends on work-hardening (measured by 7). According to 
the experiment, eqn. (2), this dependence is not stronger than linear (we 
will see below that it is found by measurement to be not less than linear 
either). Our assumption is therefore 


Tg=a+4+ fr. eee ak ahh SLAC 
From (5) and (6) it follows 


Hy — %9(G',/G) + (1 —B2)Gy/G@2— (1 —B:) ; hata) 
1—B, een 

Comparison with (2) gives «=7,(7=0) and B=dr,/d7 in terms of the 

measured quantities a and b if G(7) is known. Following Cottrell 

and Stokes (1955,) Adams and Cottrell (1955) we have defined in § 4 a flow 

stress ratio for large deformations 


Zit fmm DO aes 


To 1 
a, 641 f 
By means of (7) this turns into 
ite! any (8) 
i ele 82.0 
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Cottrell et al. also discuss a corrected flow stress ratio M = 7,6,/71G, 
which means 


1—f, - 1 — (drs/dr), 
meg Tidrcldr),. 5 Fi A . r (9) 


The deviation of M from 1 in the measurements on Al and Cu shows 
that dry/dr > 0 (see (6) ). It should be noted that the corrected flow stress 
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ratio asa function of temperature (7’,) does not directly give the temperature 
dependence of ts. From Seeger’s expression (1) it follows 
TG La ty,(d log v/dr), 


= ; 10 
7,4, 14+7,(dlogv/dr), (1) 


Me 
In a plot of M against 7’, rapid changes of 7,(7',) will show up, so that 
the corresponding critical temperatures 7’) can be determined in this 
way. However, in order to determine 7,(7’) itself it is necessary to know 
d log v/dr. 
On the basis of this interpretation of the flow stress ratio we also 
understand the temperature dependence of the rate of work-hardening 
since the relation 


+=—=1, J=— ot en ne et 
is valid for any work-hardened state in stage II (and III) where flow stresses 


and work-hardening rates are measured at two temperatures 7’, and 7’. 
To prove (11) we derive from (3) 


dr dig , dts 

da da" da 
and from (6) 

dr dtg Il 


i SAS Se 
By definition (4) is again 


drt¢ dry = 
(aa), | (Ge), = eal 


Thus with (8) one gets (11). To sum up we have shown that measurements 
of the flow stress ratio and work-hardening rates as a function of tempera- 
ture yield the same information on 7,(7') if G(7') is known}. The flow stress 
ratio is, however, much simpler to determine than the tangent of the flow 
stress-strain curves. If G(7')is not known one cannot eliminate its effect 
on the flow stress ratio with the help of 3;,(7’). On the contrary eqn. (11) 
offers a valuable cross check of the measurements by these two methods. 
In the case of our data on nickel 77/739(9y;(20°K)/9,,(7)) = constant = 1-05 
for 77=78°, 200° and 300°K{. As far as 9;,;/9;,;= constant independent 
of temperature (which is true for 78°, 200°, 300°K for nickel) the relation (11) 


+ Preliminary dynamic measurements at Stuttgart on the temperature 
dependence of the elastic moduli of some of our nickel crystals after deformation 
indicate temperature coefficients (dG/dT’)/G—=4-2 x 10-4/°K between 200°K and 
300°K and (dG/dT')/(G=2-5x 10~4/°K near 150°K. The corrected flow stress 
ratio M—=7,G1/(71@,) (fig. 19) and the corrected rapid hardening rate 9,,/@ 
(fig. 13) thus become temperature independent between about 200°K and 300°K 
while the respective rise with decreasing temperature below this range is 
diminished by 2-5 x 10-4/°K. 

t The 7,/7, data for Al (Cottrell and Stokes 1955) and Cu (Adams and Cottrell 
1955) and the scarce values on #;;(7') available in the literature (Staubwasser 
1954, Sosin and Koehler 1956, Noggle and Koehler 1957 and Berner 1957) 
indicate that (11) also holds for these metals, 


eer 
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also holds during stage III. We have, however, pointed out already in §3 
that 77/9; is abnormally high for 7'=20°K so that (11) fails here also. 
We think that the discontinuities of flow to be discussed in §6 are 
responsible for the high };;; at 20°K and may be also for a 5% too high 
a value of 9); at 20°K causing the slight deviation from (11) in stage IT 
(see fig. 13). 

The validity of the relation (11) and therefore (6) essentially shows 
that the contributions to flow stress of the dislocation cutting process, 7T,, 
and of the elastic interaction of dislocations, +,, are proportional to each 
other. This fundamental law has already been anticipated by Cottrell 
and Stokes (1955), and Seeger (1955). ; 

According to the measurements, eqn. (2) and the equations derived from 
it do not hold for stage I where presumably 7, has a major share of the flow 
stress and is not proportional to 7,,. The anomalously large change in 
flow stress for commercial nickel during stage I could mean that the 
dislocation forest is different from that of high purity nickel. There is 
also the possibility that typical alloy hardening effects come into play 
here like the necessity to destroy non-random distributions of the 
component atoms at the beginning of slip. 


5.2. Stage III 

It has been shown in § 3 that the work-hardening rate decreases at 
higher deformation as a consequence of what has been called ‘dynamical 
recovery’. The resulting stage III is approximately linear and extends 
over about 40° strain in the case of nickel. This stage sets in at a : 
critical stress 7,,;, the strong temperature dependence of which will now 
be analysed in terms of a dislocation model proposed by Diehl et al. (1955). 
These authors ascribe dynamical recovery to the escape of screw 
dislocations from the main. slip planes by means of thermally activated 
and stress-aided cross-slipt. This assumption is supported by the 
observation of cross-slip lines with the electron microscope at the 
beginning of stage III (Diehl et al. 1955, Mader 1957)i. The screw 
dislocations mobilized in this way from piled-up groups can promote 
strain and diminish the flow stress thus leading to a reduced rate of 
work-hardening. 

Figure 20 illustrates the critical step of the cross-slip process. An 
extended dislocation in the main slip plane (111), width 27), recombines 


+ In 1954, Leibfried and Haasen concluded that thermally activated cross-slip 
is responsible for the gross temperature dependence of the stress-strain curve of 
face-centred cubic crystals. The specific process considered at that time was 
the formation of slip lamellae, out of a work-hardened slip plane, at high 
temperatures and large strains. 

+ Note added in proof.—Dr. 8. Mader, Stuttgart, has investigated recently 
one of our nickel crystals of orientation No. 6 by the electron-microscopic 
techniques he has developed for copper. He found essentially no cross-slip at 
7=4-1 kg/mm? which is below our ty,=5'5 kg/mm? at 300°K. 
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over a certain length 2/, and separates again into partial dislocations 
in the cross-slip plane (111) where it is now free to move. The activation 
energy for this process has been calculated by Schoeck and Seeger (1959). 
For a stress of the magnitude of the critical shear stress 7) acting in the 
corss-slip plane the activation energy HZ, turns out to be about | ev for j 
aluminium, 10ev for copper. The difference is caused by the respective ~ 
widths of the stacking faults in dislocations in these metals which are 
determined mainly by stacking fault energies y of 200 erg/cm® for Al, 
40 erg/em? for Cu. We will try now to estimate HZ, and accordingly y for 
nickel from our measurements on stage IIT. 


eS = 


Fig. 20 
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Cross-slip of extended screw dislocation (Schoeck and Seeger 1955). 


For this purpose we make use of the stress dependence of the activation 
energy for cross-slip also calculated by Schoeck and Seeger (1955). It 
is obvious from the magnitude of H, particularly for Cu that at normal 
temperatures the applied stress has to contribute a sizeable part of the 
work to bring together the partial dislocations. The activation energy 
for cross-slip is lowered by the stress 7 according to 

E=E,— clog (7/79) o & ie © lie sa 
where c=0-35ev for Al, 3-2ev for Cu per order of magnitude change in 
stress in the interesting stress range (10-8G) (Schoeck and Seeger 1955). 
Here 7 is the applied stress times n, the number of dislocations in a 
piled-up group, i.e. 7 is the stress acting on the leading dislocation 
of the pile-up. There is considerable experimental evidence (Seeger 
1957) suggesting that n is a constant, about 20 to 30, during stages IT 
and III of work-hardening. Given the stress and the temperature the 
rate v of cross-slip is expected to follow from 

v=v,, exp(—H(r)/kT). Se es yk 

It is now assumed that the beginning of stage III is characterized by 

this rate v reaching a certain value v,;;; which is constant for different 

temperatures and metals. Then the temperature dependence of the 
corresponding stress 7;;; follows from (12) and (13) as 
kT a 


E 
log ty y.= (“Gt +loe 78) — og 2 > «(Te 
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We have already shown in §3, particularly fig. 14, that the stresses ™ 
we have found for nickel as well as those measured on copper by Blewitt 
et al. (1955) and Berner (1957) depend on the temperature in the formt+ 
logty;=A—BT. We are now able to estimate H,(Ni) and c(Ni) from the 
measured constants A and B, $3.6, when we eliminate the parameters 
log n and log (v,,/vzz;) in (14) by means of the corresponding measurements 
on Cu. It seems reasonable to assume that these parameters are about 
equal for the two metals at the beginning of their respective stages IIT. 
For 7) we have chosen the low temperature critical shear stress as a repre- 
sentative stress level at which we can compare the activation energies 
H, of the two metals (7>= 1000 g/mm? for Ni, 79=250 g/mm? for Cu 
(Blewitt), t)= 200 g/mm? for Cu (Berner)). The result of the analysis of 
the beginning of dynamical recovery, as pictured in fig. 14, in terms of 
the theory of thermally activated cross-slip, expressed by eqn. (14), is 
e(Ni) = 0-92c(Cu, Blewitt), or c(Ni) ~c(Cu, Berner) and H,(Ni) = 0-81#,(Cu, 
Blewitt) or H,(Ni)=0-82H,(Cu, Berner). 

These activation energies H, for cross-slip are (like the widths of the 
stacking faults) proportional to G'/y where G, the shear modulus, is 80°, 
larger for nickel than for copper (at 300°K). From this and from the 
measured /)’s we derive a stacking fault energy for nickel of the magnitude 

_ Ho(CujG(Ni))  _ 5. 
YNi— B,(Ni)G(Cu) You= 2"2 Yow 
Le. yy,=90 erg/em?t. This still relatively small stacking fault energy 
and the strong elastic repulsion of the partial! dislocations place nickel 
with copper in the group of metals with widely extended dislocations as 
distinct from aluminium with its narrow stacking faults§. 

There are some other theoretical and experimental published data 

concerning the stacking fault energy of nickel. In terms of electronic 


7 We have been unable to confirm the relation 7*;;;7'’=constant in the case of 
nickel. Such a temperature dependence of 7,,;; has been proposed by Friedel 
(1955) considering a breakdown of Lomer—Cottrell dislocations at high stresses. 
This mechanism did not account for the data on copper either. 

t We have tried to determine y directly from the stress 7°,;; extrapolated 
to 7’=0 where the partial dislocations must be brought together without any 
thermal activation. In the pertaining relation (Seeger 1957) 

7117 =2G(0-056 —y/Gb)/n, 
however, the second term in the brackets is small compared to the first. Thus, y 
cannot be obtained accurately from this equation (which yields reasonable 
values for 7). 

§ For Al single crystals only a few values of 7,;; have been published 
(Staubwasser 1954, Sosin and Koehler 1956, Jaoul 1957) which were measured 
at temperatures low enough to give a well-developed stage II. Evaluation of 
these data along the lines of this section yields 

c(Al) Ss; E(Al) < 0:4 

e(Cu)”. #,(Cu)~ 
Therefore y for Al is much higher than for Cu, as has been expected. The differ- 
ence between the calculated (Schoeck and Seeger 1955) and observed ratios of 
c(Al)/c(Cu) and H(Al)/H,(Cu) is attributed to the rather high values of yyy 
used at room temperature where stage II is barely existent for Al. 
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structure there is no reason why nickel, situated between cobalt 
(y= 10 erg/cm2) and copper (y=40erg/em?) in the periodic table, should 
have a high stacking fault energy, as Seeger (1955) has pointed out. The 
overlap of 3d electrons of neighboured atoms is small at the end of the 
third long period, and so is the difference in energy between the stacking 
orders ABCABC and ABABAB of close-packed planes. Actually, a 
hexagonal phase of Ni has been obtained by condensation from the vapour 
at low temperatures (Reimer 1957). Since the energy of a twin boundary 
and the stacking fault energy are closely related to each other in the 
face-centred cubic lattice, the large frequency of twins in Ni after 
recrystallization also suggests a low y. 

A more direct observation of stacking faults in nickel results from an 
analysis of shifts of Debye—Scherrer lines after filing (Paterson analysis). 


oe ae ee 


Smallman and Westmacott (1957) determine the stacking fault probability — 


of nickel to be about half of that for copper, indicating yy; 2 yc, Christian 
and Spreadborough (1956) conclude from line shifts and the change of 
electrical resistivity with deformation that nickel belongs in a class with 
copper regarding its stacking fault energy. Seeger (1955) interprets older 
data on the plastic properties of nickel to yield yy;=80 erg/em?, in good 
agreement with our estimate. 

The stacking fault energy discussed above is an important parameter 
in the theory of the critical shear stress and, in general, of the dislocation 
cutting process. This becomes more complex than indicated by (1) in 
the case of low stacking fault energy metals (Seeger 1955) and will not 
be treated here in detail. 


§6. Discontinuous Sire AND MECHANICAL TWINNING 


At large stresses, reached during deformation at low temperatures, 
Blewitt and co-workers (1957) found two types of irregularities in the 
stress-strain curves of copper single crystals. Type A, fig. 21 (a), shows 
up in the stress-strain curves, at 4:2°K, of crystals whose initial orientation 
is in the [100] corner of the stereographic triangle. Type B, fig. 21 (6), 
is observed on crystals with orientation in the [111] corner, deformed at 
4-2°K. By means of a special x-ray technique (Sherrill e¢ al. 1957) and by 
surface investigation after etching the big load drops in the B curves 
could be attributed to the formation of deformation twins, on {111} 
planes with (112) as shear direction. These twins spread across the 
specimen in the form of a Liiders band while the load fluctuates about the 
low value (2-3-4 in fig. 21(b)). The serrations in A curves could be 
shown to belong to localized avalanches of slip. At 78°K this discontinuous 
slip did not occur in crystals of any orientation while twinning has still 
been observed, near [111], starting from a neck in the specimen. In the 
meantime discontinuous slip also has been observed during low temperature 


deformation of many other polycrystalline metals (Basinski 1957, Wessel 
1957). 
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Discontinuous slip in copper single crystal at 4:2°K according to Blewitt 
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Discontinuous slip and twinning in copper single crystal at 4-2°K 
according to Blewitt et al. (1957). 


416 Peter Haasen on the Plastic Deformation 


It should be noted that Blewitt et al. (1957) found the twin plane to be 
identical with the main slip plane in most of the cases although twinning 
also has been observed on the conjugate slip plane, especially after neutron 
irradiation. By irradiation the flow stress level can be raised considerably 
so that a smaller amount of normal slip precedes twinning. 

In a note added in proof to Blewitt’s paper it is stated that Au and Agt 
have been observed to twin in a similar fashion to Cu at 4-2°K and 78°K. 
No evidence of twinning has been found, however, in the cases of Al and Pb, 
deformed at 4:2°K. Such a difference in behaviour would just have been 
expected if twins were generated out of stacking faults in dislocations, 
and if the stacking fault energy determined the stress necessary to nucleate 
a twin. Actually, the stacking fault ribbon in an extended dislocation 
in the (111) plane can be considered a mono-atomic twin layer of finite 
width. If the stress is large enough to pull away one of the partial 
dislocations against the attraction of the surface energy of the stacking 
fault the latter can spread and, by means of a spiral mechanism, grow 
into a macroscopic twin. It seemed, therefore, worth while to look for 
mechanical twinning in nickel since the magnitude of the necessary stress 
might give one more clue for the stacking fault energy of nickel relative 
to that of copper and aluminium. 

Rather unexpectedly no load drop of the kind of fig. 21 (b), indicating 
the formation of twins, could ever be observed with these nickel crystals, 
at temperatures of 4-2°, 20°, 78°, 200° and 300°kK, even for orientations 
favourable for twinning according to Blewitt e¢ a/. (1957). The stress level 
in nickel was up to a factor of two higher than in copper. Since we were 
able to reproduce twinning in Cu in some test runs, our soft machine was 
not the reason for the failure to induce macroscopic twinning in Ni. 
Furthermore, two nickel crystals were pulled in Blewitt’s hard machine 
at 4:2°K. Figure 22, Pl. 12 shows part of the chart of one of these runs 
(crystal no. 12b), indicating discontinuous flow in nickel at 4-2°K but no 
twinning load drops. Our machine did not record the small serrations but 
one could often see them directly on the pull-rod at 20°K. 

Metallographic investigation of the nickel specimens deformed at 
4-2° and at 20°K, however, frequently showed sharp bands near the 
necked down region (fig. 23 (a), Pl. 13) which were clearly different from the 
irregular low temperature slip lines present after heavy deformation. 
Polishing or etching did remove the slip lines but not the sharp bands 
(fig. 23 (b), Pl. 13). X-ray investigation, using Blewitt’s and other 
oscillating techniques with monochromatic radiation, proved much more 
difficult in Nithan in Cu due to the stronger distortion of the nickel lattice. 
Tn one case only (crystal 12 a) a reflection of a (111) twin could be identified, 
with (111), [110] as the main slip system. The twin plane in this case is 
identical with the ‘unpredicted’ slip plane (Cahn 1951) which comes 
into action at the [010]|-[011] boundary of the stereographic triangle. 
The sharp bands visible on this crystal also had the orientation of the 


} See also the recent investigation of Suzuki and Barrett (1958). 
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‘unpredicted’ slip plane while other crystals deformed at 20°K showed 
‘twin bands parallel to the conjugate slip plane (111). 

We have to conclude that in nickel it takes the stress concentration 
of a neck in the specimen to induce what appear to be twins, even at 20°K, 
while copper shows twinning clearly at 78°k. Thus, either the stacking 
fault energy y of nickel is considerably higher than that of copper, or, if 
we rely on our estimate of a rather low y from the dynamical recovery data 
($5), we are left to assume that other factors besides y also determine the 
formation of mechanical twins in face-centred cubic crystals}. We think 
that the large work-hardening capacity and the corresponding high 
dislocation density of deformed Ni impede the formation of macroscopic 
twin lamellae (see also Suzuki and Barrett 1958). Discontinuous slip 
on the other hand seems to be necessarily connected with the high stress 
level reached in low temperature deformation where at the head of only 
a few piled-up dislocations the theoretical shear strength is reached and 
new dislocations are generated freely (Seeger 1957, Haasen 1958). 
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SERVI (1958) has observed certain ridges on an otherwise smooth surface 
produced by the action of a special etching agent on zinc. He interprets 
these ridges as showing the course of dislocation lines pre-existing in the 
specimen, suggesting that accumulated solute near the dislocation lines 
has made the metal more etch resistant. No other interpretation seems 
probable, but Servi confesses to being unable to interpret certain features 
in his second picture, which appears to show a ‘ Frank—Read source ’ in 
operation. These features, he says, are the central ‘dip’ in the curve 
joining the two spiral centres, and two ‘dots’ (we should say rather 
‘dashes ’) near to it. He is in fact wrong to say that the first of these 
features is unexpected on the assumption that this is a ‘ Frank—Read 
source ’. The conditions of motion of dislocation lines in this case are 
evidently such that a forward-facing ‘ vee ’ in the line persists, instead of 
rapidly filing in (this implies strong frictional hindrance to the motion 
of these dislocation lines); ‘vees’ will then be produced at the place 
corresponding to Servi’s dip in every operation of the source except the 
first. The ‘ dots’ or ‘ dashes’, however, remain to be explained. It is 
first necessary to point out that on any interpretation the etching process 
has made gaps in lines which must in reality have been continuous. The 
best we can therefore do is to fill in these gaps in the way which gives the 
simplest acceptable interpretation of what is seen. Figure 1, traced from 
Servi’s fig. 1 (b), shows the visible lines of the pattern we have to complete. 
The first which occurred to us was that there were really two ‘ Frank—Read 
sources’ close together, as shown in fig. 2 (a). The other two fixed 
points of the sources would then be at the upper ends of the two dashes, 
which would represent the only visible portions of a short length of disloca- 
tion line in the form of a letter w. This is the interpretation which 
requires the least length of gap, but the objection to it is that each source 
should have a re-entrant ‘vee’ for itself—there should be three * dips ’ 
instead of one (fig. 2 (b)). The interpretation indicated in fig. 3 is free 
from this objection, and we can see no other unobjectional interpretation 
of comparable simplicity. In this, two dislocations emerge from each 
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centre. All the gaps are in the second member of the pair. This may be 
explained by supposing that it collected less impurity because the first 
had swept the ground ahead of it. The strong frictional drag indicated by 
the configuration of the dislocation lines is attributable to interaction with 
the solute. Then the frictional force on the second would be less than 
that on the first, explaining the unequal spacing between the two spirals. 
We think this is rather more likely than the other possibility, that the two 
dislocations are a pair of partials in a very widely extended dislocation. 
A possibility not to be excluded on the evidence is that the broken and 
continuous lines, respectively, of fig. 3, show two successive positions of the 
same dislocation lines, at each of which they have rested long enough for 
some segregation of solute to occur; so that only one dislocation line 
emerges from each centre. 


REFERENCE 
SERVI, I. 8., 1958, Phil. Mag., 3, 63. 


Evidence for the Trapping of Interstitial Atoms in 
Irradiated Tungsten 


By M. W. THompson 
Metallurgy Division, A.E.R.E., Harwell 


[Received December 16, 1957] 


THE residual resistance of tungsten specimens has been measured at 4°K 
after reactor irradiation at —196°c. Exposure in the BEPo reactor at 
Harwell (~ 101? n/em?/sec) for six days doubled the residual resistivity 
from its original value of 0-06 microohmcm. It has been reported 
elsewhere (Kinchin and Thompson 1958) that the recovery of such an 
increase occurs in two main stages of roughly equal magnitude; one 
centred at —80°c with an associated activation energy of 0-48 ev and 
one at 320°c with 1-7ev. Slight recovery, however, occurs over the 
whole range of observation from — 196 to 450°C. On the limited evidence 
of this work it was suggested that the recovery in the — 80°C stage is due to 
interstitials becoming mobile and that in the 320°c stage to the migration 
of single vacancies. Recovery in the intermediate range and at tempera- 
tures above 320°C was attributed to the release of trapped defects and to 
the migration of vacancy clusters. 

The experiment to be discussed below has provided evidence for the 
trapping of interstitials. A tungsten specimen was irradiated for | hour 
which caused an increase in residual resistance of 0-91% (+0-01%). 
Annealing periods of 30 minutes at —40, 100, 200 and 350°c respectively 
produced the recovery curve shown in the diagram. The specimen was 
heated to 1500°c which restored its resistance to the pre-irradiation value, 
it was then irradiated for 12 hours which caused a resistance increase of 
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881°). Annealing was carried out as before and this recovery curve is also 
shown in the diagram. The most obvious difference between these results 
and those reported previously is the absence of any large change in 
resistance after annealing at —40°c. It is also clear that although 
increasing the irradiation time to 12 hours causes more recovery at 
100°o and 200°¢ there is still very little change at —40°c. 

These results may be explained by postulating the existence of several 
types of trap to which interstitial atoms could become attached (Lomer 
and Cottrell 1955). Such traps could presumably take the form of 
substitutional impurity atoms whose size is less than that of the tungsten 
atom. <A rough lower limit to the concentration of traps may therefore 
be taken as the concentration of impurity atoms; this quantity may be 
deduced from the spectroscopic analysis and the residual resistivity to be 
about 10-4. Now after a 1 hour irradiation the observed increase in 
resistivity indicates a concentration of displaced atoms of the order L0~*. 
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(Ry —R,)/(R;—Ry) versus annealing temperature. Ry,—residual resistances 
after annealing at temperature 7’. R o—=residual resistance before 
irradiation. A;—residual resistance immediately after irradiation. 
Dashed lines indicate the position of principal recovery stages observed 
after 6 day irradiations. 


It is clear therefore that there are insufficient interstitials to fill all the 
traps. Under these conditions mobile interstitials will have a far greater 
chance of becoming trapped than of annihilating vacancies or of migrating 
out of the lattice. The slight recovery observed at low temperatures may 
be understood by supposing that interstitials released from post-irradiation 
positions by thermal excitation are immobilized again by the plentiful 
supply of traps and that such a process only causes small resistivity changes. 
The recovery observed at 100 and 200°c after the 12 hour irradiation may 
be explained by assuming that in this case there are insufficient deep traps 
to hold all the interstitials and that when migrating at these temperatures 
they will start to annihilate vacancies or will migrate out of the lattice, 
either process giving rise to larger changes in resistance, 
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The following question now arises: is the previous conclusion that 
interstitials are immobile in the lattice until —80°c correct, or are they 
mobile and become trapped at the irradiation temperature? In either 
ease the mechanism of changing traps, proposed above, will explain the 
experimental results provided that the ‘frozen in’ interstitial is regarded 
as being a special type of trapped interstitial. However, if the former 
answer is correct then the amount of recovery at —40°c should reach a 
limiting value for long irradiations which produce sufficient interstitials to 
occupy all the traps. On the other hand if the interstitials are immobile 
at —196°c then their number should continue to increase until the 
concentration of defects is high enough to cause instantaneous recovery 
of further damage. For practical irradiation times in BEPO the recovery 
at —40°c would not show saturation in this case. Further experiments 
are proceeding in which recovery after longer irradiations will be studied. 
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REVIEWS OF BOOKS 


Ferroelectricity in Crystals. By Haren D. Macaw. (London: Methuen and 
Co. Ltd.) [Pp. xi+220.] 27s. 6d. 


FERROELECTRICITY is a branch of physics that has developed very rapidly in — 
recent years. Until 1935 only one ferroelectric was known ; now many are — 


known, and more are being discovered at an increasing rate. The volume of 
literature is so large that it is hard to keep up with it, and anyone starting work 
on the subject must find it very difficult to get an adequate picture of its present 
state. A general review of the subject was therefore badly needed, and it 
could hardly have been better done than in this book. It is surprising how 
much information has been put into 200 small pages. Little of importance 
seems to have been omitted, though naturally for a fuller development of some 
subjects the reader must refer to the original papers. This conciseness, without 
loss of clarity, has been obtained by careful arrangement and the author’s sound 
judgment of what are the most significant facts in published work. As might 
be expected from Dr. Megaw’s own interests emphasis has been placed on the 
crystal structure of the substances described, but this is quite right since an 
understanding of the nature of ferroelectrics must ultimately depend on a 
detailed knowledge of their atomic structure. Naturally the greater part of the 
book is devoted to experimental work, since in this subject theory has lagged 
far behind experiment ; but three chapters of the book are devoted to theories 
of ferroelectricity, with an assessment of their merits. 

It is to be hoped that this book will have the success it deserves, and that it 
will be possible to bring out second and subsequent editions at not too long 
intervals. Ay ees 


Constantes Sélectionnées : Diamagnétisme et Paramagnétisme, Relaxation Para- 
magnétique. By G. Forx, C. J. Gorter and L. J. Smits. (Paris: Masson 
& Cie.) [Pp. 325.] 8,800 fr. (clothbound) ; 9,700 fr. (stiff covers). 

Apart from a brief section on paramagnetic relaxation (by Gorter and Smits), 

the tables deal mainly with susceptibility and provide a nice example of the 

well known distinction between the physicist, who makes accurate measure- 
ments on materials of unknown purity, and the chemist, who makes measure- 
ments of unknown accuracy of pure materials. The lack of cooperation between 
physicists and chemists is illustrated by the great rarity of entries in heavy type 

(indicating great reliability) or asterisks (indicating agreement between two or 

more authors to better than 1°). Probably the main value of tables such as 

these lies not so much in the numerical data they provide but in the biblio- 
graphy ; here the bibliography is admirably organized, but is somewhat in- 
complete in the field most familiar to the reviewer. Dice 


Viscous Fluid Flow. [1—Turbulent Flow. By Sura-I Par. (Princeton, U.S.A.: 
D. Van Nostrand Co., Inc.) [Pp. 277]. 48s. 

Tuts book is an attempt to provide a survey of what is known about turbulent 
flow, presumably for American engineering students of the subject at post- 
graduate level. It is comprehensive so far as titles of chapters are concerned, 
but the treatment of the different aspects of turbulence is rather scrappy and 
superficial. The author has a bias towards those parts of the subject that lend 
themselves to mathematical description, and this leads him into making a 
curious and misleading selection of material. The more physical aspects are 
presented poorly, or left out altogether, This is not a book that the present 
reviewer would recommend. CK as 
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J.T. BARTLETT and J. W. MITCHELL Phil. Mag. Ser. 8, Vol. 3, Pl. 9. 


Platelets of gold produced by exposing and heating a silver chloride crystal 
containing 0-05 mol, aurous chloride. (x 875.) 


Fig. 2 


Particles of gold with associated systems of decorated prismatic dislocations 
in a silver chloride crystal. The photomicrograph is taken in a (111) 
plane. (x 875.) 
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Hexagonal network within a silver Dislocations interacting with imper- 
chloride crystal. (x 875.) fections in a recrystallized silver 


chloride crystal. (x 875.) 


The dislocations of tilt sub-boundaries and three dimensional networks 
within a well annealed silver chloride crystal. (> 875.) 
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Fig. 5 


Dislocations introduced by straining a recrystallized silver chloride crystal. 
(x 875.) 


Fig. 6 
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Dislocations in a silver chloride crystal strained after recrystallization. 
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Fig. 22 


Elongation ~ 


j 
| 


) 


Discontinuous slip in nickel single crystal 12 at 4:2°K. 


PETER HAASEN Phil. Mag. Ser. 8, Vol. 3, Pl. 13. 


Fig. 23 (a) 


Nickel crystal 22c after 54% strain at 20°K. Magn. x 8-5. 


fickel crystal 22c¢ after 54°, strain at 20°K, polishing, and etching. Magn. x 8-9. 


Neck at the left side of the picture. 


